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Servo-Stabilization of Low-Frequency Oscillations 
si rey in Liquid Propellant Rocket Motors 


By FRANK E. MArsie, Pasadena, California!) 


Introduction 


Two points of view may be taken with regard to the undesirable pressure 
oscillations in rocket motors which arise from instability of combustion, ac- 
_ coustic resonance, coupled oscillations of chamber pressure and propellant flow 
- rate, as well as from more obscure sources. One is to eliminate the underlying 
cause of instability through change in mechanical design or modification of 
_ propellant properties; this is possible when the mechanism of instability is 
understood and its removal is not detremental to rocket performance. The 
alternative is, as was demonstrated by H. S. TsIEn [1] 2) to modify the system 
_ dynamics by utilizing a feedback servocontrol which, for example, senses press- 
_ ure fluctuations in the combustion chamber and modifies the propellant feed- 
ing rate at the proper frequency and phase to damp the fluctuation. Servo- 
- stabilization provides the distinct advantage that stability need not be a major 
_ concern during rocket design, relying upon the feedback system to insure stable 
operation. Furthermore these concepts suggest the possibility of eliminating 
~ empirically an undesirable oscillation even when its basic cause is not known. 
The purpose of the present paper is to demonstrate the process of servo- 
stabilization by illustrating how the requirements may be met in the specific 
~ case of low frequency oscillations in liquid propellant rocket motors. The advan- 
tage of this example is that there exists a well developed theory for these 
oscillations, as shown, for example, in the discussions of M. SUMMERFIELD [2] 
y and of L. Crocco [3], and hence the removal of instability may be demonstrated 
. analytically. The analysis of low frequency instability in monopropellant and 
bipropellant rocket motors is reviewed employing the Satche diagram intro- 
duced by TsIEN [1] for stability analysis of systems with time lag. Then it is 
shown in detail how utilizing the same diagram, the response of a simple feed- 
back servo loop may be prescribed so as to assure stability for all values of the 
_ propellant time lag. It must be kept in mind, however, that the specific nature 
of the example places no restriction upon the generality of the principle. 


+ 
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1) Daniel and Florence Guggenheim Jet Propulsion Center, California Institute of Technology. 
2) Numbers in brackets refer to References, page 34. 
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It has been the author’s privilege to prepare this extension of Professor : 
H. S. Tsren’s original work on servostabilization of rocket motor with the 


benefit of frequent discussions and suggestions from Professor TsIEN. 


_. Dynamics of the Liquid Propellant Rocket Motor 


In reviewing the principles of instability in liquid propellant rocket motors 
it is well to consider immediately the bipropellant motor since the simple 
monopropellant rocket appears as a special case, the same stability criteria 
applying after suitable simplification. Because only low frequency oscillations 
are under consideration, the gas pressure at any time will be assumed uniform 
over the rocket combustion chamber. If the dependent variables of the system 


are taken to be the chamber pressure, the oxidizer flow rate and the fuel flow _ 


rate, three equations relate these quantities to each other and to the physical — 


q 


configuration. They are the continuity relation for the combustion chamber and 


two conditions of dynamic equilibrium for the fluid in the oxidizer and fuel 


lines. Following Crocco’s notation [2] let M, be the mass of burned gas in the_ 
combustion chamber, m, the mass rate of burning, and m, the mass rate of 


ejection from the discharge nozzle. Then the relation 


ap (Ma) = My —'m 


(a) 


expresses the conservation of mass. Suppose the state of the dynamic system is | 
always close to an equilibrium value defined by M,, m, = m, = m, then the | 


residence time which, on the average, the burned gas spends in the chamber, 


ise M,|m. Anticipating the use of perturbation analysis, the fractional — 


variation of burning rate and discharge rate may be introduced, 


y= = and ee 


Me — m 
m m 


together with a dimensionless time z = ¢/0, to write the mass conservation as 
a ({M,\ + 
Bes a) 


The propellant burning rate m, differs from the total propellant injection 
rate m, because of the combustion time lag. The time required for the transfor- 


mation of propellant into products of combustion may be represented by means 
of the integral 
t 


/ te. T,) at’ = C = const (2) 


t-—t 


———— 
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where /(p, T,) is a function depending upon the mechanics of the process. The 


_ chief quantities upon which this function depends are #, the uniform chamber 
_ pressure and 7,, the ambient gas temperature in the neighborhood of the 
injector. The meaning of this integral may be clarified somewhat by thinking 


of / to be, for example, the rate of heat transfer to the liquid propellant, a cer- 
tain total quantity C of heat being required for the complete transformation of 
liquid propellant to vapor. The principal quantities associated with such a heat 
transfer process are the ambient temperature and pressure, so that the time t 


_ required for transformation will be large or small depending upon whether the 
_ heat transfer rate is small or large. Now if the time lag were independent of 


time, the fractional variation of burning rate uw, would be identical with the 


_ fractional variation of injection rate measured at the time t earlier. Denoting the 


dimensionless time lag t/0, = 6, constant time lag would imply j,(z) = u,(z — 0). 
However, if the time lag is increasing during the process, the burning rate is 


_ depressed below the injection rate in proportion to the product of the propellant 


flow rate and the rate of change of time lag. Therefore using dimensionless 
quantities and neglecting a term of second order in the small variations, 


pio(2) & wi(z— 8) — $2. 3) 


_ The dependence of the time lag upon time is, in reality, governed by equation (2) 
and differentiating this with respect to time gives, 


fle, TO) — #166 — 2). Ta] {1 — 2] =0. ( 


Since it is assumed that the chamber pressure and ambient temperature vary 


_ only slightly from their steady state values, # and T,, the values of / may be 
- approximated, following Tsren [1], through a Taylor expansion about = ?, 


- 


; 1, = fT, 
flo), Tol] 1B, T) + MP tp —F] + MP [0 ~T, 
flo —2), e— o] 1@, T) +") p0— 9-2] 


f Upon substituting these relations into equation (4) and deleting second order 
differences, the variation of time lag may be expressed 


4 dx _ dlogf(p, Ty) le sald Mes a lass l 


at dlogp 


P 


Tf) (Gt-») - Be 5 
ee jo 


dlog T, 


g 
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For a monopropellant rocket the combustion gas temperature is nearly constant, 
so that the second term on the right side of equation (5) vanishes. Calling 
dlogf(P, T,)/Ologp =, TsteN’s generalization of Crocco’s treatment is ob- 


tained. Considering the physical processes involved, it is likely that the value ~ 


of n is less than 0-85. For the bipropellant rocket the combustion gas tempera- 
ture may experience considerable variation, resulting from the variation of 
propellant mixture ratio, and consequently another parameter 


dlogf(P, T,) 
dlog T, 


=u 


enters the problem. Denoting the fractional variation in chamber pressure as 
gp =(p— )/p and the fractional variation of injector end temperature as 
y= (1,= T,)/T,, the expression for dimensionless time lag variation is 


do 


© ws nl plz — 0) — g@)) + ave — 9)-= v@)]. (6) 


The fractional variation in total propellant injection rate u;, which occurs | 


in equation (3), is most conveniently expressed in terms of the individual 
fractional variations of oxidizer flow, (uo = (Mp) — M)/mp, and of fuel flow, 


he =n; = my) lm, and the steady state mixture ratio 7 = Mo|My. Again follow- | 
ing Crocco in defining H = (7 — 1)/[2 (ry + 1)] depending on the steady state 


‘mixture ratio, it is an elementary calculation to show that, to the first order, 


it 1 
Mi= (5. ar H) [lo + (> = H) Ly. 
Consequently the propellant burning rate is just 


ole) = (3 + H) wole— 8) + (7-H) we 0) | 


— n[p(z — 6) — (2)] — gly — 4) — y@)]. 


In calculating the variation of mass flow m, from the nozzle, it is convenient 
to consider the nozzle flow as quasi-steady, that is, the flow corresponds to a 
time sequence of steady equilibrium flow patterns. Then the ratio of instanta- 
neous to steady mass discharge rate is 
m p [( Ty \1/2 
= eRe (=*) ) (8) 


ge 


(7) 


where 7, is the gas temperature at the exhaust end of the combustion chamber. 
TsIEN has shown [4] that the quasi-steady approximation is adequate for the 
low frequency oscillations under discussion here. For the monopropellant 


‘ 


ch: 
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rocket, where 7,, is invariable, the mass discharge rate depends directly upon 
the combustion chamber pressure # and the fractional variation of mass dis- 
_ charge rate is simply , = y. Where the variation of mixture in a bipropellant 
rocket causes appreciable variation of temperature this too must, according to 
equation (8), be accounted for in computing the rate of gas discharge. Suppose 
the gas temperature at the exhaust end is 7,, at the time ¢. Recalling that the 
gas temperature depends upon only the propellant mixture ratio, the temperature 
T,, can be related to the injection rates of propellant which formed the gas 
volume. The generation of gas took place at a time 6, earlier then the gas 
appeared at the nozzle end of the chamber, that is at a time ¢ — 0,. Further- 
more the time required between injection of the propellants and generation of 
the gas is the time lag r. Therefore the gas temperature at the exhaust end of 
the chamber is that corresponding to the mixture fatio at the time ¢—1— 0, 
or at the dimensionless time z — 6 — 1. Now, if 7 is the instantaneous mixture 
_ ratio, the change of gas temperature from its equilibrium value is, approximately 


a lr 7. 
- Furthermore the variation of mixture ratio from its equilibrium value follows 
- directly from the definition 
r= —° wr (1+ Mo — Hr) 
Denoting the dimensionless slope of the gas temperature curve, dlog T,/dlogy, 
_ by 2K, the effect of combustion gas temperature upon fractional variation of 


mass discharge is K [ju(z — 6 — 1) — wr(z — 6 — 1). Thus the complete ex- 
pression is 


fe(2) = 92) — K[polz — 8 — 1) — wz — 6 — 1). (9) 


From this calculation of the effect of mixture variation upon gas temperature 
it is clear that y= 2K (uM — my), and consequently equation (7) may be 
simplified by writing the term 


ple — 8) — ple) = 2K {[jolz — 6) — myle — 8)] — [eol2) — a2)]}- 


Finally the variation of gas storage within the combustion chamber is 

~ easily found by computing the variation of M, = af o, dV with time. Since the 
gas density @, is directly related to the temperature and uniform pressure in 
the combustion chamber through the equation of state, the ratio of the actual 
- mass M, to that value M, which would be there under steady conditions may 


_ be written be is 
M, Mee dV _1y yee dV 
SS te = aa 4 Be eS 
uM, ape, V [ p(2)] Peay 
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Neglecting the product of small variations in g(z) and the integral, 
a (M,\  , 4 (i av 
da Gz) Tee da Ve 3h Al en 
Now the variation of the integral arises simply from the fact that gas of temper- 
ature T,, is being ejected from the exhaust end of the chamber while gas at 
temperature T,, is being generated at the injector end. Thus the variation of 
the integral is simply the difference in the integrand at the two ends multiplied 
by the mean volume flow rate 
af Ee W (Bo 
edie Re Pe edi 


gt ge 


The gas temperatures T,, and T,, depend upon the local mixture ratios; more 
precisely, T,; corresponds to the mixture ratio injected at the time z — 6 while 
T,. corresponds to that injected at z—6—1. Then utilizing the previous 
calculation of the variation of gas temperature with propellant mixture ratio, 


it follows directly that 
d (M,\_ dg | 
dz (iz) as 


(10) 
PAK MG on ee 8 Tae eee 


Introducing equations (7), (9), and (10) in the original expression for the mass 
conservation in the combustion chamber, equation (1), gives this relation in 
terms of the desired dependent variables 


#4 (Ln) pt nge—9) 


0. 
(11) 


This description of rocket dynamics may be completed by relating the 
individual propellant flow rates to the proportions of the propellant supply 
system and the properties of the propellants. Consider the supply for each propel- 
lant to consist, Figure 1, of feeding mechanism (e.g. turbopump), a propellant 
supply line, an injector nozzle, and a variable capacity situated in the line so 
close to the injector nozzle that the inertia of fluid between the capacity and 
the injector nozzle may be neglected. Since the process is the same in oxidizer 
and fuel lines, notation differentiating the two lines will be postponed. Denote 
by subscript 1 conditions in the propellant line just downstream of the pump 


+{ Kwe-0~1) + [2K@-1)— } A] ple 9) ~ 2K 4 n9(2)| 
z [3K E20 1) [2K@-1) + = —H| ple 3) — 2K 4 u,(2)| = 


wh 


b 
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and by 2 conditions just ahead of the injector. Following TsIEN [1], the variation 
_ of m,, the pump discharge rate, from its equilibrium value may be written 
peat Se eae) (12) 
my : Pr 
When % =oo this relation indicates no response of mass flow to a change of 
line pressure, that is, a positive displacement system. Conversely when « = 0, 


OXIDIZER PUMP 


ROCKET CHAMBER 


FUEL PUMP 


Figure 1 


Schematic diagram for general bipropellant rocket. 


no change of line pressure may be effected by changing the flow rate, which is 
- characteristic of the constant pressure feed or pressurized tank. Intermediate 
values of « correspond to particular types of pump and if the response time of 
the turbopump system is very long with respect to the period of chamber 
_ pressure oscillation, the pump will operate at constant speed and « corresponds 
- to the slope of the pump characteristic (Figure 2). The circumstance where the 
rotative speed of the pump may not be considered constant has been treated 
poy Y.C..Les, M.R. GorE, and C.C. Ross [6], in which another degree of 
- freedom appears in the problem, although this feature has little influence on 
_ the result. 

The pressure difference ~, — p along the propellant lines is due primarily 
to the liquid accelleration through the injector nozzle and the inertial reaction 
of the fluid column in the propellant line. The frictional loss will be neglected. 
If S is the area of the propellant jet issuing from the injector and g the propel- 
lant density, the pressure difference , — p across the injector m?/(2 @ S?). 
- Furthermore if/and A are the length and cross sectional area of the connecting 
line, the inertial pressure difference follows from the momentum equation 


(61 — Bb.) A= (Umm). 
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The mass flow m, through the injector may differ from the mass flow m, in the 
connecting line because of the variable capacitance situated between the two. 
Since the propellant is considered incompressible, the difference of the two 


values of mass flow may be expressed m— m, = dC/dt where C is equal to the | 


instantaneous fluid mass contained in the capacitance. This capacitance may 


be represented conveniently in the dimensionless form K = C/ (m 6,). Then if — 


PUMP DISCHARGE PRESSURE, p, 


CONSTANT SPEED 
PUMP CHARACTERISTICS 


m, 
PUMP MASS FLOW RATE, Mm, 


Figure 2 


Approximation to constant speed pump characteristic. 


the pressure drop across the propellant injector under steady operation, 
m?|(2 @ S?), is denoted 4, two dimensionless parameters may be introduced; 


the inertia of the connecting line J =/m/(2 Af A 6,), and the injector 
pressure ratio P = /(2 Ap). Finally eliminating variable conditions at the 
pump and at the capacitance, the dynamic equilibrium of a propellant line may 


be written 
du 


2 
T+ C41 wt Pets 


x = 
ace tS 


ax 
rata (13) 
where (P + 1/2) « has been denoted &. An equation of this type holds indepen- 
dently for the oxidizer line and the fuel line. 

In [7], R. SABERSKY has discussed the behavior of a simple system where 
the propellants may not be considered incompressible and consequently wave 
propagation along the propellant line must be accounted for. The result is, of 
course, that the effective inertial constant J of the propellant lines decreases 
as the number of nodes in the propellant line increases. 
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Stability of the Monopropellant Rocket 


| When the rocket utilizes only a single propellant and the propellant lines 
_ are of fixed capacity, the system is described by two independent variables ¢(z) 
and y(z) representing the fractional variations of chamber pressure and fuel 
flow. The dynamics of the monopropellant rocket follows from the general 
‘description by allowing the mixture ratio to vanish so that H > —1/2, K +0, 
and one of the variations in propellant flow vanishes. Then the two differential 
- equations are, for the combustion chamber 


da 
“a + Lm) p +m ole —8) — we — 6) = 0 hd 
and for the fuel line 
, d = 
Pao ee +l) pe PoH 0. eu 


"The stability or instability of the system, that is whether the amplitude of 
particular oscillating solutions will diminish or grow, is most easily discussed 
_ through introducing the Laplace transform, defined as 


[o,2) 


q 3 0(s) = | 9 e-** da (16) 


0 


_ for the pressure fluctuation g(z). Similarly defining M(s) as the Laplace trans- 
_ form of y(z), each of the differential equations may be transformed, as shown 
in the Appendix, into inhomogeneous algebraic equations in the complex 
- variable s. The inhomogeneous portion depends upon the initial conditions of 
the problem and hence is associated with transient solutions. The corresponding 

_ homogeneous equations 


(s+1—n+ne**) G(s) —e °° M(s) = 0, | 
POs) +(JstC+1I)M)=0 | 


- describe the free oscillations which are of interest in the stability problem. For 
‘nontrivial solutions the determinant of the coefficients of equations (17) must 
vanish, that is 


(17) 


|stil—n+ne ES gr? 
| =G'(s)=0. (18) 
| P eae Gedy | 


Equation (18) is a transcendental equation in the complex variable s, the roots 
of which determine the stability of rocket chamber pressure oscillations: If any 
root of G’(s) = 0 posesses a positive real part the system is unstable. To check 
‘on the existence of such roots, it is possible to apply the familiar Nyquist 


:| 
| 
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criterion which traces the course of the complex function G’(s) as s describes a 
counterclockwise contour in the right half of the s plane consisting of the 
imaginary axis and a large semicircle in the right half plane. Then the negative 
change in argument of G’(s) is the difference in number of zeros and poles of 
in the right half s plane. Graphically this number appears as the number of 
complete revolutions made about the origin by the trace of G’(s). If the number 


of poles is known independently, the number of zeros in the right half s plane | 


is known and the stability of the system is determined. 


The present problem is complicated by the time lag 6 which not only makes | 


| 
| 
| 
| 
| 
| 


calculation of G’(s) laborious but accentuates the fact that the value of 6 is | 


usually not known accurately. To improve this situation TsIEn [1] has intro- 


duced another technique based upon a suggestion of SATCHE [8] which sepa- 


rates the function G’(s) into two parts, the first containing e~ °° 


as a factor, the — 


second independent of the time lag. In the present example the separation may — 
be effected by noting that the equation (18) may be expressed as the com- — 


bination of two determinants 


istl—n 0 | 
D,(s) = | (19) 
| n —] 
D,s) = | (20) 
| P fS Got | . 
such that equation (18) may be written 
D,(s) + e~*° D,(s) = 0. (21) | 


Following TstENn in calling g,(s) = e~‘° and defining g,(s) = —D,(s) /D,(s), the 
characteristic equation becomes 


81(s) —-g,(s) = G(s) = 0, (22) 


where the complex function G(s) is represented as the difference of the two 
vectors g,(s) and g,(s). 


p 


Assume for the moment that G(s) has no poles in the right halfplane. Then 


if s describes its previous contour the number of zeros which G(s) has in the 


right half-plane is equal to the number of complete clockwise rotations made — 


by a vector having its head on the trace of g,(s) and its tail on the trace of 
inp : te 

£o(s). Since all values of g,(s) lie on or within the unit circle, the vector G(s) can 

make a complete rotation if the trace of 2(s) either passes through the unit 


circle or encircles the origin, Conversely if the trace of g,(s) lies outside the unit _ 


circle and does not encircle the origin, the system is definitely stable. Actually 
the function G(s) may possess poles in the right half-plane due to zeros of the 
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determinant D,(s) which constitutes the denominator of g,(s). Thus the number 
of clockwise turns which g,(s) makes about the origin must be modified by 
-adding the number of poles induced through the roots of D,(s). This number is 
found by applying the Nyquist criterion to D,(s) directly. 
To illustrate the technique consider the injection rate to be constant so that 
M(s) = 0; instability under these circumstances has been denoted intrinsic 
instability by Crocco [3]. The determinants degenerate; D,(s) =s+1—n 
while D,(s) = 1, so that 


gs) = -(~=") - =. (23) 


The curve g,(s) reproduces the contour of s in the left half-plane with a scale 
factor 1/n and displaced a distance — [(1 — n)/n] to the left of the origin. 
Clearly g,(s) misses the unit circle when (1 — )/n > 1, that is forO Sn< 1/2; 
the Satche diagram is shown in Figure 3 for m = 1/4. When > 1/2 the possi- 


Figure 3 


Satche diagram illustrating intrinsic stability, ™ = 1/4. 


12 FRANK E. MARBLE _ ZAMP | 


bility of instability shows itself through the intersection of the two curves, 
Figure 4, so that clearly » = 1/2 is the critical value of the index n describing 
the pressure dependence of the propellant time lag. This criterion for possible 


intrinsic instability may, due to its simplicity, be obtained analytically with — 


Figure 4 


Satche diagram illustrating intrinsic instability, n = 3/4. 


almost equal ease as the graphical solution. However the advantages of the 
technique employing the Satche diagram become more apparent when applied 
to more involved problems. As a second example consider the combined oscil- 
lations of chamber pressure and propellant injection rate which occur when 


Gel kPa Sy hee 


DN) w 


Then the general expressions for D,(s) and D,(s), given by equations (19) and 


: 


: 
| 
: 
: 
: 
; 
: 
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(20), lead to the relations 


3 D,(t w) = 3 (1 — n) —4074+[4 (1 —n) +4 3Jio 
and 
Dio) =3(n+ 5) +[4n]io. 


The curve of g(t w), Shown in Figure 5 for 2 = 0-2, indicates no intersection 
with the unit circle. To close the Satche diagram, take s = R e’® where 6 de- 
‘creases from 2/2 to —z/2. Therefore for large values of R, go(s) behaves as 


Figure 5 Figure 6 
Satche diagram for stable monopropellant Satche diagram for unstable monopropellant 


tocket; J = 4:0, P = 3/2,a=—1,% 0-2. rocket; J = 4:0, P = 3/2,a = 1, 7 = 0°6. 
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so that the g,(s) curve closes by a large arc in the left half of the Satche diagram, 
progressing clockwise. Since g.(s) neither crosses the unit circle nor encircles 
the origin outside the unit circle, instability can occur only if a zero of D,(s) 
exists which cancels a root of g,(s) in the right half s-plane. Since D,(s) has a 
zero at only an imaginary value of s, this is not possible, and consequently the. 
system is certainly stable. 
When the propellant time lag index m is increased to a value of 0-6, the. 
situation is considerably modified for, as shown in Figure 6, the possibility of : 
low frequency oscillations is indicated by the intersection of g,(7 m) with the. 
unit circle. The Satche diagram closes in precisely the same manner as for the 
lower value of m. Therefore a rocket and propellant having parameters with 
these values will exhibit a low frequency instability for some range of values _ 
for the propellant time lag. The method for determining the critical time lag, 
above which the instability will occur, has been shown by TsIEN [1], in his. 
original discussion of this graphical analysis. | 


Stability of the Bipropellant Rocket 


The behavior of the bipropellant rocket motor is complicated by the addi- 
tional degree of freedom associated with the second propellant line and by the 
periodic changes in combustion gas temperature caused by variations in mixture 
ratio. The system is described by three independent variables, (2), Mo(z), and | 
J4;(z) representing the fractional variations in chamber pressure, oxidizer and 
fuel flow. The three differential equations are, for the combustion chamber 


d 
Poy 
dz 


(1 —n) p+ (z—) | 

ae [3K pale — 8 —1) + [2K @—1)— 4 —H] po(z — 0) | 
a 2K 4 ule) u. (3K mee 0-1) + [2K @-1)+4 —H] | “ 

X pi(z — 6) — 2K a 1(2)\ =0 | 


and for the oxidizer and fuel lines, 


di Bs a 
Jovae + CotN mot he=0, fy F+4+ (G+) m+ Be=0, (25) 


Introducing the Laplace transform as before, defining M,(s) and M,(s) as the : 
Laplace transforms of Holz) and «u;(z), each of the differential equations may be 
transformed as shown in the Appendix. The corresponding set of homogeneous 


> 
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algebraic equations is 
[st+1—n+ne °*|D(s) 


4 ae [Ke*+ 2K q-1)— 4 —H]—2K a} MIs) 


(26) 


| —{e*? [Ke'+2K@-1)+ 5 —H] —2K q| M,(s)=0, 


Py P(s) + Jos + Co +1) M,(s)= 0, 

P, B(s) + (Jy s + C4 +1) M,(s) = 0., 
_and the stability of the system is determined from the roots of the characteristic 
equation constructed from these coefficients. Defining two third order deter- 


minants D,(s) and D,(s) as the logical extension of those given in equation (19) 
and (20), 


sti-n 0 =2Kq 
D,(s) = Po Jost Cg+1 2Kq , (27) 
P, 0 Ipstbp41) 


| n Ke‘+(2Kq@-1)-H- 5) —[Ke‘+ (2K (q@-1)+H-5)] 


Ex Jos +o +1 0 hae 


P, 0 Jes + Cy +1 


The characteristic equation is identically in the form given by equation (21), 
_and the stability analysis follows in an identical fashion. 

To discuss the stability of any bipropellant rocket two sets of parameters 
must be known; H, K, n, and g which depend upon the propellants and their 
mixture ratio, and Pj, P;, Jo, Jy, %», % Which are associated with the propellant 
pumps, feed lines, and injector. In particular the parameters H, K, and probab- 
ly » and gq are not independent of each other; H depends upon the propellant 
mixture ratio and K upon the slope of the combustion temperature versus 
mixture ratio curve at the steady mixture ratio under consideration. For the 
examples to follow, these values will be taken as given in Table I for various 
mixture ratios. 


Aa nT) re be 


Table I 
7 | a (TR | aT ,|dv | H K 
2-75 5065 0 0.2335 | 0 
2:50 5010 360 0-2125 | 0.09 
2-25 4880 640 0:1925 | 0-1475 


2-00 4706 850 0-167 0-189 
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For the first examples choose the mechanical properties of the system such that 


Oe Opa IOs Po= Peels TREE 1 | ee 


and consider cases corresponding to two values of the steady state mixture © 


ratio and two values of the coefficient 1 = dlogf/dlogp, the value of q is 
considered to be zero. Then the values of the functions D,(¢ w) and D,(t w) may 
be computed according to the relations 


| D,(i @) = {= (Lm) Ee Taek Ces n)| ot 


< (p25Pigs 4 
+iwl + Z (1 — n) oe ke 
Dio) =|F n+ 5 +5 Kosin —3n ow 


ie ie (eee eae K cosa}. 
The curves of g,(z w) for n = 0-2, calculated using these relations, are shown 
in Figure 7 for 7 = 2-5 and in Figure 8 for 7 = 2-0. They are nearly identical 
despite the considerable difference in mixture ratio; neither of them intersect 
the unit circle. Proceeding as for the monopropellant rocket, it follows that for 
large values of R, g,(Re'°) behaves as — R e'®/n so that the curve closes by 
a large arc in the left half of the Satche diagram, progressing clockwise. Further- 
more the associated Nyquist diagram for D,(s), shown in Figure 9 for n = 0-2 
and vy = 2:50, behaves as n J; J) R® e?"® for large values of R. Therefore this 
diagram closes with a large arc of 2 a radians starting from the upper branch in 
Figure 9 and progressing clockwise to the corresponding lower branch (not 
shown). Since the origin is not enclosed, the system is certainly stable. The 
same conclusion may be reached for the lower mixture ratio, 7 = 2-0, n = 0-2. 
The stability curve for the corresponding example with 1 = 0-6, shown in 
Figure 10 for 7 = 2:5 intersects the unit circle. Therefore a rocket propellant 
combination having parameters of these values will exhibit a low frequency 
instability for some range of time lag 6. As in the first two examples, the closing 
arcs in the Satche diagram and the associated Nyquist diagram for D,(s) indicate 
no further instability of the system. 

Because the mechanical properties of the two propellant feed lines are 
nearly identical in the foregoing examples, the change in mixture ratio asso- 
ciated with an oscillation in chamber pressure is very small and the behavior is 
much like that of a monopropellant rocket without line elasticity. Consequently 
the Satche diagrams are particularly simple. The influence of mixture ratio 
oscillation becomes important when the two lines have different inertial con- 
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‘ 

- 
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Figure 7 


Satche diagram for stable bipropellant rocket; 7 = 2°5, 1 = 0°2, Gp = &, = 1:0, Po = Py = 10, 
Jo ec 2-0, J+ = 15. 


stants and injection pressures. As an example consider a rocket for which 
my = oy == 1-0; joe #0; jy= 1-0, P= 10, P= 4-0, 7=2:0,n=0-2. 
The equations for D,(¢ w) and D,( w) may be written 
‘ 55 49 al 
D,(i o) = Fy ee [+4 0—m) 0 


49 


+iofS 42-404}, 


ZAMP VI/2 
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Figure 8 
Satche diagram for stable bipropellant rocket; 7 = 2-0, n = 0:2, % = &-= 1:0, Py = Py 
Fo= 2:0, Jp = 18, 


= 1:0 


? 


Depa [33 n—4o%n+ 29 ex_meK (S cosw + 15a sina} 
: 49 Wt = !) : 
+to 2 ht = IS (2 ye (15 COS@ =~ sino). 


The Satche diagram for this example, F 


igure 11, exhibits large loops that orig- 
inate in the trigonometric terms of D, 


(? w). These terms arise because the change 
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Figure 9 
Nyquist diagram for D,(s) associated with Satche diagram of Figure 7. 


in combustion temperature associated with a change in mixture ratio does not 
affect the chamber discharge rate until the residence time 6, later; consequently 
this disturbance has a reduced time lag of unity. The physical phenomenon 
underlying the loops is easy to see. When a gas mass of low temperature, caused 
by a fluctuation in mixture ratio, flows to the chamber discharge nozzle the 
mass flow rate increases [c. f. equation (8)] causing the chamber pressure to 
drop. This reduction of chamber pressure induces a change in propellant flow 
rate with a consequent variation of mixture ratio which may lead to another 
negative fluctuation in gas temperature. At certain frequencies these. oscilla- 
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Figure 10 
Satche diagram for unstable bipropellant rocket motor; 7 = 2°5, 1 = 0-6 Xp = %, = 1:0 


Po = Py = 1-0, Jg = 2-0, Jy = 15. 


tions will reinforce each other so that they will, as in the case of the present 
calculation, be slowly damped. As pointed out by Tsren, the Satche diagram 
indicates possible slowly damped oscillations when the £2(2 w)-curve passes 
close to the unit circle. Furthermore it is clear that there will be an infinite 
number of frequencies at which such a reinforcement will take place. In the 
present example only two show up to any extent. The amplitude of an oscilla- 
tion which originates in this manner depends to a considerable extent upon the 
response of the propellant feed system. If the change in pressure produces a 
large change in propellant flow, the oscillations may be severe; if the response 
is poor, the oscillation will disappear quickly. Now at high frequencies the 


7 eet aes eee ee 
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Figure 11 


Satche diagram for stable bipropellant rocket motor; 7 = 2-0, n = 0°2, %& = & = 1-0, Py = 1-0, 


P, = 40, Jo = 4:0, Jz = 1-0. 


inertia of the propellant supply causes the amplitude of its response to decrease 
and consequently this particular resonance disappears. Mathematically this 
may be seen in the fact that for large values of w the term Jo J+ @? governs 
the behavior of D,(s) whereas the trigonometric terms are at most proportional 
to o. 

The details of the Satche diagram near the unit circle, Figure 12, show that 
the curve lies outside the unit circle. The closing arc of the g,(s)-curve is the 
same as those discussed in previous examples. The Nyquist diagram of the 
denominator D,(s) is shown in Figure 13, indicating no zeros of this determinant 
in the right half of the s-plane and, as a consequence of this and the above 
considerations, the system is conclusively stable for all values of the com- 


bustion time lag. 


Be. Frank E, MARBLE ZAMP 


Figure 12 


Details near the unit circle for Satche diagram of Figure 11. 


Servo-Stabilization of Low-Frequency Oscillations 


Where possible unstable operation of a rocket is indicated, stability may 
often be assured by modifying the mechanical design of the propellant supply 
system or injector, or in the extreme, by change in the chemical propellants. 
Alternatively the stability may be assured by introducing a feedback loop 
connecting the combustion chamber and the variable capitances in each pro- 
pellant line. In the simplest case for the monopropellant rocket, Figure 14, the 
loop will consist of a pressure measuring device attached to the rocket chamber, 
an amplifier, and an appropriate servomechanism to actuate the capacitances 
in the propellant lines. With this scheme the capitance x will vary with time 
in a manner prescribed according to the variation in chamber pressure. Employ- 


ing a dimensionless form, the correspondence between capacitance and chamber 
pressure variations may be written 


x(2) = F( 7) 92), (29) 
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Figure 13 
Nyquist diagram for D,(s) associated with Satche diagram of Figure 11. 
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where F(d/dz) is a linear differential-integral operator. This relation simply 
indicates that the capitance and the chamber pressure are connected through 


SERVO 
CHAMBER 
PRESSURE 
PROPELLANT be 


Figure 14 


Schematic diagram of monopropellant rocket incorporating feedback system. 
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a linear differential relation with constant coefficients, the form of which de- 
pends upon the pressure pickup, the amplifier and the servomechanism. The 
problem is then to design the feedback circuit, that is to determine the form of 
F(d/dz), such that stable rocket operation is achieved. In a practical case this 
resolves to the design of an amplifier to work with a given pressure pickup and 
a given servocontrol so as to produce the required form of F. 

An adequate form for F can be determined only through a stability analysis 
of the entire rocket system with feedback loop, the behavior of this system is 
described by equation (14) for the rocket chamber, equation (13) for the pro- 
pellant line, and a relation of the type equation (29) for the feedback loop. The | 
homogeneous portions of the Laplace transform of these equations are 


[stl—n+ne**] g(s)—e °° M(s) =0, | 
P ols) +[Js+o+1] M(s) +[Js?+ 6s] K(s) =0, | (30) 
F(s) p(s) — K(s) =9, | 


where K(s) is the Laplace transform of x(z), the variable capitance; F(s) is the 
transfer function of the feedback loop. The stability of the system is dependent 
upon the roots of the third order determinant constructed from the coefficients 
of equations (30), 


Sell arena aes? 0 | 
P Jst+C+1 Sfs*+Es1. (31) 
| F(s) 0 —1 | 


| 


It is most convenient to expand this in terms of second-order minors using the 
last row for the expansion. Then defining a new second-order determinant 


| u 0 . 
D(s) = | : | (32) 
Js+O41 Js?+s | 


the characteristic equation obtained by setting the stability determinant 
[equation (31)] to zero is just 


D,(s) + e-*°|D,(s) + F(s) Dis) |=05 
where D,(s) and D,(s) are those corresponding determinants which are appro- 
priate for the monopropellant rocket without the feedback loop. Then calling 


oO — ache! D,(s) c 
8e(s) = Di) + Fi Dis” (33) 


the stability analysis is reduced to that corresponding to equation (22). 
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Stable rocket operation is to be assured by choosing a form of the feedback 
transfer function F(s) so that the Satche diagram satisfies the criteria previously 
discussed. The direct procedure discussed by TsIEN [1] is to choose a simple 
stable Satche diagram and solve for the feedback transfer function. In some 
cases, in particular for the bipropellant rocket, this procedute leads to a form 
for F(s) which is not a rational function and hence is not conveniently reali- 


_ zable. The basic difficulty here is that an excessive amount of information has 


been prescribed by fixing the entire Satche diagram, the necessary conditions 
for constructing a stable Satche diagram are actually simpler. Clearly it is 
necessary to be assured that g,(0) lie to the left of the unit circle, that is 
&2(0) < —1. Then since the curve g,(7 w) has always a vertical tangent at w= 0, 
it is necessary that the curvature be sufficiently small that g,(i w) does not 
pass into the unit circle as s increases. In most cases, this may be guaranteed by 
leaving the curvature and scale of the Satche diagram at w = 0 unchanged by 
the feedback circuit. Finally it must be assured that the diagram does not 
misbehave for large values of w; it must not loop about so that it passes through 
or around the unit circle at high frequencies. This condition may be satisfied in 
general by arranging that the feedback system leaves the original Satche 
diagram uninfluenced for large frequencies. This is easily done by designing the 
amplifier to be unresponsive at high frequencies, that is the amplifier shall ‘cut 
out’ at high frequencies. These features may be incorporated into a rather 
simple transfer function for the feedback loop which is sufficient to eliminate 
instability. 

From equation (32) it is observed that D(s) behaves as s for small values of 
s so that in order to change the value of g.(0) by choosing F(s) it is necessary, 
according to equation (33), that F(s)~1/s for small values of s. Furthermore 
the conditions that the scale and curvature of the Satche diagram be unchanged 
near w = 0 requires that in the power series expansion for F(s) D(s), the coef- 
ficients of s and s? vanish. These conditions may be met by choosing F’(s) to be 


of the form 


1 a+bs (34) 


eS « i 0% 2)? 


F(s) 
where a, b,c are constants and F(w)~1/m* for large values of m and hence 
cuts out with sufficient rapidity. Then clearly the expansion for 


a+06s) (Js+ ¢) 
(1 Ss c2 $2)4 


D(s) F(s) = 4 


can have the coefficient of s vanish only if a + bs is proportional to the conju- 
gate of Js + ¢, that isa+bs=f(—Js+ ¢) where / is a new constant. Then 
furthermore the expansion for 


D(s) Fs) =f Gaye 
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can have a zero coefficient of s? only if —J? + 2c? ¢?=0 so that the constant 
c is just 
paaretalie (35) 
y2e 
Finally the number / may be fixed by setting the value of g,(0), for at s = 0 


DON, es (1 — ) (¢ +1) 
£2(0) = D0) + fC? ow (C41) — PHO’ 


Thus if it is desired to make g,(0) a certain prescribed value —y where y > 1, 
then the final constant f is just 


Bey Saw tae EG 1 ee ' 
} ; aaiesea Go) 


For the unstable monopropellant rocket considered, where «= 1, P = 3/2, 
J = 4-0, and = 0-6, the actual value of the feedback transfer function i 
easily found to be ; 
25 1 — 2:03 s? 
LSet spe eee 


where y has been arbitrarily set at the value 1-50. The corresponding Satche 
diagram, Figure 15, indicates that the rocket is stable against low-frequency 
oscillations for all values of the time lag. 

For the bipropellant rocket the general approach to the problem of servo- 
stabilization is the same and the actual mechanisms are only slightly more 
involved. In the most general case, two variable capacitances, Figure 16, are 
involved and it will be assumed that each is provided with its own feedback 
transfer function. 


nol) = F(-5) of). (37) 
“g(@) = E(4-) ot) (38) 


and the problem is then to determine forms of fy and F; to assure stable oper- 
ation of the rocket by submitting the full set of equations describing the 
chamber pressure, dynamic equilibrium of each propellant line, and each feed- 
back circuit to a stability analysis. As shown in the Appendix, the set of alge- 
braic relations obtained after application of the Laplace transform is 


Bee 9 pet ot ee ee | 
—e*|_3Ke '— (-2K+H—5)] Mjs) =0, | Ce 
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Figure 15 


Stable Satche diagrams for.monopropellant rocket obtained using simple criteria on feedback loop. 
Py y(s) + [Jos aes Col M,(s) + [Jo s?+ lp s] K,(s) =0, 
P, os) + [Ips +1 + S,] Mis) + [Fp 8? + & 8] Kils) = 0, 
Fx(s) p(s) — Ko(s) =0, 
F,(s) y(s) — K;(s) = 0, 


(39) 


where K,(s) and K,(s) are the Laplace transforms of the variable capacitances ; 
F,(s) and F;,(s) are the transfer functions of the feedback loops. The stability 
of the system is determined by the roots of the fifth order determinant construc- 
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Figure 16 


Schematic diagram of bipropellant rocket motor incorporating a general feedback system. 


ted from the coefficients of equations (39), 


stl—nine s® ¢86 [Ke*—(2K+H+5)] | 
ty Jos+1+hq 
P, 0 | 
EAS) 0 
F,(s) 0 
Ls poate : 1 (40) 
et [Kent — (2K+H- >)| 0 0 | 
0 Jops?+igs 0 | 
Jps +1 4%, 0 st +Gs\- 
0 —1 0) ) 
0 0 of | 


This determinant is most conveniently expressed through expansion by second- 
order minors, employing the last two rows for the expansion. Defining two new 


: 
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Bsird-order determinants 


wv 


4 ss 1 

Ke (2kK+H+ ) 

E eo ee eet et 

ne 

; ‘ (41) 
; a [K e-*— (2K+H-)| 0 

.. 

4 

; 0 Yi ee 

, 

{ Jp7st1+, 0 

Sx. 1 

at (2K+H+ 5) 

i ES) Fistt+y | 

0 

. (42). 
— |Ke-s— (2k+H-=)| 0 

. 0 ys? +hys 

str, 0 

: 


_ the characteristic equation obtained by setting the stability determinant [equa- 
_ tion (40)] to zero is 


D,(s) + e~*° [D2(s) + Fo(s) Dols) + Fs) Dy(s)] = 0. (43) 


- The transfer functions of the two feedback loops, /(s) and 1;(s) appear as 
- multiplicative factors; each of the determinants depends upon only the original 
_ rocket configuration and is independent of the feedback circuit. The new form 
_ of the function g,(s) is 


—D,(s) Fe) 
Das) Fils) Dols) + Fils) Dyls) ’ na 


so that clearly as much flexibility may be expected by stabilizing with one 

propellant line as with both. Therefore it will be supposed, arbitrarily, that the 

- feedback servo will vary only the flow of oxidizer. This means that F(s) = 0 

in the present formulation so that the problem becomes very similar to the 
- monopropellant case and may be treated in the same manner. 

From the expression for D,(s), equation (41), it appears still proper to 

choose F(s) to behave as 1/s for small values of s and hence it is appropriate 
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to assume again fi 
Ks) = > qe | (45) 


ce [1 —‘c? s?]? 


together will the criteria that for small values of s, the scale and curvature of 
&o(s) remain unchanged. It is sufficient then to make the coefficients of s and s?_ 
vanish in the expansion of [%(s) D,(s) for small values of s. For small s, the’ 
determinant D,(s) may be expressed in the series 


D,(s) = d s+ d2 24d) s84-.., (46), 
where 
dp =(K+H+ >) oo (1+q), | 
af) =(K+H+ 5) [h+) + pol+ Koo +), | en 
a) =(K+H+ 5) hh+K[hO+G) + bol-7RGQO+G), 


ne “ . torecte z 
dy =Khhk-7K[h (1 + oy) + Col te Koo 0 + G,) : 


Then as for the monopropellant rocket, the dependence upon s may be elimi- 
nated by writing a + bs = f (d‘) — d® s). Then F,(s) D(s) becomes 


Fis) Dos) = 1 (2 = a9) (at + a? s+ ap s+ ++) 


(1 — c? s?)2 s (48) 

so that in order for the coefficient of s? to vanish, c2 must be chosen | 
oA ya! 2 ae 

c= 2 bea, aa ? (49) 


Then, finally, the value of # must be chosen so as to insure that 
LOH =amt 


and this is determined by equation (44), to be 


so that 
i 1 
f= — gaye [Pol0) —— Dy). (50) 
These results then determine completely the form of the feedback transfer 


function required to stabilize a liquid bipropellant rocket against low-frequency 
oscillation of the type described, 
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Consider, as an example,’ the unstable bipropellant rocket with parameters 
having the values, 
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Cer = pare p= 20, Jp lo5) f=—29; Danes 


which was analyzed ascend and whose patcne ete is shown in here 10. 
_ From these values it is easily determined that 


a =301, d?=6-15, d) — 2.89 
and ‘consequently | . ? 
eZ a=301, b=-6-15, c?=1-60. 
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Figure 17 


 Satche diagram for stabilized bipropellant rocket obtained using simple criteria for feedback loop. 
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Finally the value of { may be calculated from equation (5) 

f= —0-62, 
so that the final form for the transfer function of the feedback loop is 


0-62 (3-01 — 6-15 s) 
Fy tt 1:6 354 


Fy(s) = (51) 


The modified Satche diagram may then be computed from equation (44) using 
F,(s) =0 and Fj(s) as given above. The resulting diagram, Figure 17, shows 
the rocket to be definitely stable against oscillation of the type discussed. 


Concluding Remarks 


The analysis which has been presented demonstrates that a quite simple 
feedback servo loop may be devised to insure stability for any case of low 
frequency or ‘chugging’ oscillation in monopropellant or bipropellant rocket 
motors. Likewise it may safely be inferred that a similar solution may be 
achieved in other cases of coupled oscillations, even though the design of an 
appropriate feedback circuit may have to be approached empirically. In fact 
it must be emphasized that the general principle of servostabilization remains 
regardless of the detailed mechanism of instability, so that even if the present 
theories of low frequency instability must be modified in the light of further 
experience, servostabilization of chamber pressure oscillations is still appro- 
priate. At frequencies above the ‘chugging’ range, when wave propagation 
along the length of the rocket chamber must be accounted for, the problem of 
sensing the. oscillations will require more care than that which has been de- 
scribed. Likewise the control of propellant injection may become a problem when 
the frequencies are increased. However, it is not strictly necessary to actuate the 
variable capacitance at the same frequency as the oscillation to be stabilized, 
for the oscillation may be suppressed effectively by damping it only every 
two or three cycles. Thus by extending the concepts of H.S.Tsren to such a 
nonlinear feedback loop, the frequency response of the capacitance may be 
relaxed considerably. Furthermore the fact that values for the pressure sen- 
sitivity of the propellant time lag are not generally known is no obstacle, 
since the feedback circuit may be designed quite conservatively without diffi- 
culty. It appears, therefore, that servostabilization must be considered as a 
natural component of the rocket, amalgamated with the control system. 
Realization of this potential benefit depends essentially upon proper engineer- 
ing balance between mechanical simplicity and electronic elegance. 
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Appendix 
Laplace Transform of the Rocket Stability Equations 


The equation describing the combustion chamber pressure oscillation is 


9 [K nol — 6-1) -(2K+H+ 5) mole 9) (52) 
= [K ule - ete yee (2K+H- >) sz — 0)] =9, 


_ while the equations of motion for the liquid in the oxidizer and fuel lines are 


d ae aed. 

Io Po 4 (Cotl) wot hots eC, a U,, (53) 
dz dz az 
d the d 

itt bet i+ Both ge th Ge =9- (54) 


Denoting the Laplace transforms of (2), Mo(2), [,(2), %o(2), and x,(z) by D(s), 
M,(s), M,(s), Ko(s), and K;(s) respectively, the transformation of equations 
(52), (53), and (54) is straightforward except for the terms involving a time 
lag. In general the transform of f(z — 0) will be 


[te — 8) e8* dz et f(z — 6) e-S@-9 d(z — 6) 
0 —0d 
ae (sy ten tf fea, 


—6 
where F(s) is the transform of /(2). As is well known, when the variable is 


retarded the transform is modified through multiplication by an exponential 
factor e~ °° and by the addition of an integral extending over the period of the 


time lag. The equations then transform to 
1 
[stl—n+ne**] Ds) + pate [« Out (2K+H+5)| M,(s) | 


(55) 
Sain ff lake (2K+H— 5) M,(s) | 


ZAMP VI1/3 
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0 0 
= 9(0) —1” e* f g(é) e 8° dé —Ke0+ [Tuglé) — py(E)] e** dé 
5 


-1-6 


+ (2K +H) e-** [ Tuglé) — mye] e-*? af (55) 
= 


0 
+ Fe" [ Luglé) + wylé)]e-* ab, 
6 


[Jo(s) +1 + £5] Ma(s) + Py ®(s) + [Jo s? + by s] Kols) 


(56) 
= Jo Ho(0) + Jo-—G2(0) + Los + £0] ¥0(0) 


[h(s) +1 + 6] MG) + BO) + ps8 4+ Gs] K; 


(57) 
= Jy us(0) + Jp 2-0) + Lys + Gs] 40), 


where ¢, and ¢, have been employed to denote 


[Ps ++ al % and [B -F 5| Xe 
respectively. Stability analysis requires only the left-hand sides of these equa- 
tions. Transient behavior of the rocket system, however, is determined by the 
complete equations where the right hand sides, in usual problems, consist in 
the initial conditions of the dependent variables. For systems with time lag, 
such as the liquid propellant rocket, equations (55), (56), and (57) show that 
in addition to initial conditions evaluated at z= 0, certain features of the 
behavior of the system must be known for a period preceding z = 0 equal in 
duration to the time lag. Physically this is clear, for due to the combustion 
time lag the variation of chamber pressure g(z) can not be described for z> 0 
unless the behavior of the propellant flow is known for the period— 6d Xz50, 
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Zusammenfassung 


Die Theorie der kleinen Frequenz-Instabilitat von Raketenmotoren mit fliis- 
sigen Treibmitteln wird iiberarbeitet, wobei im besonderen die Methode der 
Satche-Diagramme erértert wird, welche durch H. S. Tsien fiir Stabilitatsrech- 
nungen linearer Systeme mit zeitlicher Verzégerung eingeftihrt worden sind. Das 
Prinzip der Stabilisierung durch Benutzung einer Gegenkopplungsschleife wird 
besprochen, und es wird gezeigt, dass schon eine sehr einfache Lésung des Gegen- 
kopplungskreises die Stabilitat von fliissigen Einzel- und Doppeltreibmitteln in 
Raketenmotoren versichern kann, unabhangig von der Grésse der Zeitverzége- 
rung. Auf Grund des beschriebenen speziellen analytischen Beispiels wird gefol- 
gert, dass solche Stabilisierungsmassnahmen fiir eine viel grdssere Klasse von 
Instabilitaten in Raketenmotoren als die betrachtete empirisch verwirklicht 
werden k6énnen. 
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Die Rolle des Adiabatenexponenten 
bei chemisch veranderlichen Gasgemischen 


Von IRENE SANGER-BREDT, Paris?) 


Ubersicht 


Der in der klassischen Gasdynamik so fruchtbare Begriff des «Adiabaten- 
exponenten» als des Verhiltnisses der spezifischen Warmen bei konstantem 
Druck und bei konstantem Volumen verliert bei thermodynamischen Str6- 
mungen bisweilen an Sinn, da die in den Strémungsgleichungen vorkommenden 
einheitlichen x-Werte sich bei Wegfallen bestimmter einschrankender Voraus- 
setzungen als durchaus verschiedenen Ursprungs und verschiedener Grdsse 
erweisen, wortiber man sich bei der Handhabung dieser sehr gelaufigen Formeln 
nicht immer Rechenschaft gibt. 

Unter «einschrankenden Voraussetzungen» sind hier die Bedingungen kon- 
stanter spezifischer Warmen der einzelnen Gase und konstanter Partialdruckver- 


haltnisse der Gasgemische verstanden. 


1) Arsénal de l’Aéronautique, ChAatillon-sous-Bagneux (France). 
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Solange sich die Stromungstechniker vorzugsweise mit kalter oder schwach 
aufgeheizter Luft befassten, waren beide Einschrankungen berechtigt und sehr 
zweckmiassig. Bei hohen Strémungs-Mach-Zahlen der Luft oder bei str6menden 
Feuergasen jedoch, besonders in den Expansionsdiisen von Raketen oder in 
durch Kernenergie aufgeheizten Maschinen, trifft keine der beiden Vorausset- 
zungen mehr zu. Die durch Vernachlassigung dieser Tatsache verursachten Ab- 
weichungen der Rechnung von den wirklichen Verhaltnissen liegen bei Tempe- 
raturen iiber 2000°K, insbesondere bei kleinen Gasdrucken, bereits ausserhalb 
der Grenzen des Zulassigen. 

Von solchen Gesichtspunkten ausgehend, ver6ffentlichte G. KLOBE {1]*) vor 
wenigen Jahren fiir iibliche Feuergase mit Temperaturen bis zu 3500°K eine — 
im Strémungsabschnitt zwischen Ofen und engstem Diisenquerschnitt global 
brauchbare — Naherungsmethode, wie sie zur Behandlung von Strémungsfra- 


gen in Raketen nach Art der A4 verwendet wurde. Bei dieser Darstellung war _ 


unter anderem die Veradnderlichkeit der spezifischen Warmen vernachlassigt. 

In der folgenden Abhandlung wird nun versucht, den Einfluss jeder der 
beiden einschrankenden Zusatzbedingungen systematisch und getrennt zu er- 
fassen sowie anschliessend verallgemeinerte Gleichungen aufzustellen, die auch 
bei Wegfall einer oder beider Einschrankungen giiltig sind und aus denen die 
klassischen Strémungsgesetze als Sonderfalle hervorgehen. Untersucht werden 
wiederum nur Strémungen idealer Gase. Es werden fiir differentielle Entspan- 
nungsgefalle streng giiltige Isentropenexponenten abgeleitet, deren sachgemasse 
Beniitzung die formelle Weiterverwendung der klassischen Adiabatengleichun- 
gen (und damit aller daraus abgeleiteten Strémungsbeziehungen) auch bei 
chemisch veranderlichen Strémungen gestattet. Es wird ferner gezeigt, dass 
diese Isentropenexponenten nur im Grenzfall unveranderlicher Partialdruck- 
verhaltnisse mit dem Verhaltnis x = C,/C, tibereinstimmen. Endlich wird noch 
die Anwendbarkeit klassischer oder halbklassischer Adiabatenexponenten in 
Fallen bestimmter, nicht im Gleichgewicht befindlicher Strémungen erlautert. 

Die Ergebnisse der Untersuchungen werden an dem einfachen Beispiel eines 
nach dem Reaktionsschema X,=@ 2X veranderlichen Gasgemisches mit den 
Zahlenwerten des Wasserstoffs veranschaulicht. 


I. Vollig eingefrorene oder «klassische » Strémung 


In ihren Beziehungen fiir isentropische Entspannung setzt die klassische 
Strémungslehre bekanntlich temperaturunabhangige spezifische Warmen G3 
bzw. C, und unveranderliche Partialdruckverhaltnisse p,/p oder Molbriiche n,; 
der strémenden Gase voraus. 

Unter diesen einschrankenden Zusatzbedingungen gibt es in der allgemeinen 
Gaszustandsgleichung drei verdnderliche Zustandsgréssen: 1. den Druck 


) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 65. 
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_ p{kg/m2], 2. die volumspezifische Masse oder Dichte [kg s?/m*] und 3. die Tem- 


ee ee 


ee ee ea Yee ae 


peratur J [° K]. 


Jeweils zwei davon, in zweckentsprechender Wahl, kénnen als unabhiangige 
Veranderliche angesehen werden. Ausgehend von der Zustandsgleichung eines 
Gemisches von z idealen Gasen, 


eer 


ee [kg/m*], (1) 


bzw. von einer einzelnen Gemischkomponente mit dem Partialdruck ~,, dem 
Molekulargewicht M; und der partiellen') Dichte 9,, 


lore at 


Pi 2 M, a2 (1 a) 
worin 
& Pi yy 
Pi mh 
o= See M= ye M; und C= 03 ?; ap 
eae: 
Gesamtdruck, mittleres Molekulargewicht und Dichte der Gasmischung bedeu- 
ten, lautet das vollstandige Differential: 
dp dT.) de , 
p Tint, ie (2) 
bzw. 
dp, a7 | da 2a 
Esl gi Ory - 


Ausserdem wird bei konstantem Partialdruckverhaltnis die Differenz der 
molaren spezifischen Gemischwarmen: 


z r by 
it (/ C,,4T — f C,,47)| 
Pee ok Slew 


4 
d(H—U) i=1 ‘ 
al = af (3) 


OG, = 


= eg SG Cam Sh AR= AR [kcal/°kmol]. 


R [mkg /°kmol] bezeichnet hierin die allgemeine Gaskonstante, 
H [kcal/kmol] die Enthalpie, 

U [kcal/kmol] die innere Energie und 

A [kcal/mkg] das mechanische Warmeaquivalent. 


1) Das heisst die Dichte, die das Gas hatte, wenn es im eingenommenen Raum v = M/g@ 


allein ware. 


| 
” 
j 


Die Entropie eines Gasgemisches in einem offenen Kessel bei konstant blei- _ 
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bendem Kesseldruck # betragt: 
?, T,e p,T,o be m4 Se 
ee oe een) | 
0 0 ne Din; M, 2 DM M; | 
o=1,, i=1 
worin | 
P ites tap 
S,= Se + 0, [fF —-AR|S ed | 
0 0 . 
: or es A | 
Srvc, Ean] a frearime tH) 
. 0 o 
ist}). 


Hierin bedeutet der Index 0 den « Basis»-Zustand bei #; = 1 at, T = 298,2°K | 
und 9; = M,/24,79 kg s?/m*; worin M; das Gewicht der Gaskomponente 7 in 
[kg/kmol] bezeichnet. Unter 


Som $0. A RIn( SE) 


ist die Integrationskonstante verstanden, das heisst der konventionelle Entro- 
piewert fiir ideales Gas im « Basis»-Zustand. 

Die Bedingungen, unter denen die Entropie konstant bleibt, ergeben sich 
durch Nullsetzen des totalen Differentials von Gleichung (4), also 


ti Le rhe an a (05e 0S; 
dS» 7p = Sn M, [( oT i di ( a ap| 
i=1 ‘ : t (5a) 
ayy, i 
2 [4.0 Par 2] -0 


1) Die molare Entropie S ; jedes einzelnen Gemischbestandteiles 7 wiirde — bei gleichen Zu- 
standsgrossen ~, T und @ wie as des puanenk? — fiir sich allein 


* fee @ Nie do” 
S$; =S0, +S, ] = > Pant SGas cof ar | 2 
p zt a Gk Q 
0 0 


betragen. Gleichung (4a, b) meee Me bei Vermischung unter Konstanthaltung des Ge- 
z 


samtvolumens vor und nach der Mischung, » % =v, auftretende zusatzliche Gemischentropie 


z t=1 
As* = —AR > % Inn,. 
t=1 
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oder 
PT 
Wee eee OS; aT 4 OS; Fi 
Ty oT us e 
L1G 
i=1 (5b) 
ny; 
= _[+¢,,4-4R-2]=0 
1: My; : 
| a 
_ oder 
a MechieoS as 
de ee (PA et (252) ae] 
7 nM; 
: t=1 (5c) 
Dy My 
E betas) oa Serco, A liei0, 
7 ny M pe ie 
7 am : 
3 
; jenachdem man (f, T), (9, T) oder (4, @) als unabhingige Veranderliche wahlt. 
_ Die Bedingung konstanter Gemischentropie ds = 0 enthalt ohne weiteres die 
. Bedingung adiabatischer Zustandsinderung dQ = dU + A ada = 0, welche — 
 allein betrachtet — nicht unbedingt reversible Prozesse vorschreibt. In den 
 Gleichungen (5a) bis (5c) betrug die differentielle Zunahme der inneren Energie 
: pro Kilomol des Systems: dU = C, dT; die differentielle Enthalpiezunahme: 
_ dH=C,4T und die bei einer differentiellen Zustandsinderung pro Kilomol 
 geleistete reversible Arbeit dA = p dv = — (RT/M) do/e; dQ ist die dem System 


pro Kilomol zugefiihrte Warme. 
Aus den Gleichungen (5a) bis (5c) folgen die Beziehungen: 


a dp _ Cy aT 
: -~ AR "Ge (6a) 
| don Cy dT 
? ff Gk OF (6b) 
und 
dp Cy dg 
sass ee 


j z a 
wenn man die Ausdriicke C, = 3m; C,, und C, = pa n; Cy, als mittlere spezi- 
i=1 i= 
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fische Warmen des Gemisches pro Kilomol einfiihrt. Diese drei Differential- 
gleichungen lassen sich im Fall konstanter spezifischer Warmen integrieren, so 
dass man unter Verwendung von Gleichung (3): 


bT “= const, (7a) 
o TU») — const (7b) 

und 
)o ”=const (7c) 


daraus erhalt. 
Dabei hat man fiir das Verhaltnis der spezifischen Warmen 
(OH/OT), T(dS/0T), C 


OUl0T) an i (OSL), C 


p 


v 


die Abktirzung x eingefiihrt. Man pflegt x, da es die Beziehungen zweier Zu- 
standsgréssen bei adiabatischer Zustandsaénderung bestimmt, als den Adiaba- 
tenexponenten zu bezeichnen. Alle adiabatisch isentropen Zustandsaénderungen 
zwischen # und T, @ und T sowie zwischen # und @ sind reine Funktionen von 
x, wenn sowohl spezifische Warmen als auch Partialdruckverhiltnisse in der 
Gasmischung fiir konstant angesehen werden kénnen. 

Diese Art von Strémung diirfte unter Gleichgewichtsbedingungen nur bei 
tiefen Temperaturen, bei denen die Veranderlichkeit der spezifischen Warmen 
vernachlassigbar wird, und bei gleichzeitig hinreichend hohen Drucken, bei 
denen die Dissoziation verschwindet, verwirklichbar sein. Sie ist auch dann 
vorhanden, wenn eine Entspannung so schnell vor sich geht, dass weder der 
veradnderliche Anteil der spezifischen Warmen noch die Dissoziation den Zu- 
standsdnderungen zu folgen vermégen. In diesem Fall sind die im Kessel herr- 
schenden spezifischen Warmen fiir die Entspannung verantwortlich. Die Stré- 
mung ist jedoch unter diesen Umstanden zwar noch adiabatisch, jedoch nur 
scheinbar isentropisch, da die Differentialgleichungen (5a) bis (5c) dann nur 
noch ein Rechenverfahren darstellen, das eine Zuordnung verschiedener Zu- 
standsgréssen bei verschiedenen Entspannungsstadien erméglicht, die wirk- 
lichen Entropiedifferentiale aber entsprechend den spiiteren Ausfiihrungen des 
Abschnitts III gebildet werden miissten. 


II. Str6mung mit eingefrorenen chemischen Gleichgewichten und 
veranderlichen, tragheitslosen spezifischen Warmen 


Wenn nur die erste der beiden eingangs genannten einschrankenden Zusatz- 
bedingungen hinfallig wird, das heisst, wenn man temperaturabhangige spezifi- 
sche Warmen mit standiger, unverziiglicher Gleichgewichtsanpassung annimmt, 
so kann man unter Bentitzung der Planckschen Zustandssumme Z bekanntlich 
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[2], [3] die molaren freien Energien F, freien Enthalpien G, Enthalpien H, inneren 
Energien U, Entropien S und damit auch die spezifischen Warmen C, und C, 
der einzelnen Gemischkomponenten 7 als Funktionen ihrer Zustandssumme 
-ausdriicken: 


B= (FE) = ART InZ;; (8a) 
G— (CS )gon = -ART (InZ;—1) = AR + ART; (8b) 
O(InZ, 
U, — (US oe = +ART? 1, (9a) 
Hy — HE") = +ART{T[“Ge?] +1{=AU,+ ART; (9b) 
4 0(T InZ, 
S;— (SP")ox = +4R [S24], | 
(10) 
: Mazi? 6 ow, 
=AR{inZ;+F [Sr] =p tp | 
sowie endlich 
: d(InZ,) O{T [d(InZ,)/OT]o} ] | 
oy ART || vaehaal oT la us 
und 
d(In Z,) O{T[O(InZ,)/0 T]o} Mi 
eS arr {Por ] + =< + i} Com PARs. Ma) 


In diesen Gleichungen beziehen sich die Konstanten (FS*)oon, (GP )ox, 
= (HS) ox und (S$*)o-x auf gasférmigen Zustand bei T= 0°K; sie 
enthalten also unter anderem die gesamten Umwandlungsenergien vom festen 
auf den gasférmigen Zustand. 
Die Bestimmungsgleichungen (6a), (6b) und (6c) fiir isentrope Zustands- 
anderung kénnen mit Hilfe der Gleichungen (11) und (12) dann in folgender 
Form geschrieben werden: 
; 
; 
: 
‘ 
: 


ag = Sin,{{1+ 7 (2m) + [eet loan Zoo Tie | oe yee} 


EE 


; pine, O{T [0(InZ,)/0T]o} ]\ 47 
ee he | oT al i 


|, a [er (eta ea/oTe | | | is 


Pair tihe | | ea Raat ta o8 h(t) 
p a, Tf] eee alt ie} | 


oT 
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Beispielsweise fiir ein einheitliches, zweiatomiges Gas mit vollangeregter Ro- 


tation, einer einzigen Normalschwingung und mit vernachlassigbarem Beitrag 
der Elektronenanregung zur spezifischen Warme, dessen Zustandssumme also 


SRTr [ ,Os/T @p\-1 D oc AS a RT \5l2 e / 
Z=222=|7—\((4%) ies) So = E | 

e Os/T Or\~} 22 RT\3/2( M foe . 

= |r| [(4) [Gea (Ar) Q | 

ware, lauteten im besonderen die einzelnen Entropieanteile der Schwinguae=s 
Rotations- und Translationsfreiheitsgrade : | 


2s/ @,/T : 
§ = AR {ins tet) (16a) 
R 
S=AR{-Ingk + 1| (16b) 
und | 
Tr (2 « M/g)? (e R)52 5 
S,= AR [ins + ZIT — Ing} (16¢) 
in Druckabhingigkeit bzw. 
tr (2a RA (e Mig 3 
5, = ARI NE + ZnT —Inq| (16a) 


in Abhangigkeit von der Dichte dargestellt. 
Damit schrieben sich die Gleichungen (13) bis (15) in diesem Sonderfall: 


‘ s 
dp [7 0,\2 e Os/T ae 7 ieee Oh dS 
pee E + Leal ee ge apa aw (3a 
Ss 
Agim Ps @,\2  e %IT PP Ree pakas 
nc: [a+ ped eet — 7 | To ee eae (4a) 
und 
ess Oe eee 7 av, as 
GP TORE Sn a ar ae eae 
POE ONE SOO ele Si eee a oiat oa 
2 owes ) (e IT — 4)2 fata PRES 


In der Zustandssumme bedeuteten: 
M das Molekulargewicht des Gases, 
h das Plancksche Wirkungsquantum, 
N, die Zahl der Gasteilchen pro Kilomol, 
e die Basis der natiirlichen Logarithmen und 
Op bzw. Os die charakteristischen Temperaturen fiir die 
Rotations- bzw. Schwingungsanregung. 
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Die Gleichungen (13a) bis (15a) sind integrabel: 


s 
Inf —In(T") + a = const, (13b) 
$ 
Ing — In(T*?) + 7, = const (14b) 
und 
s 
7/5 2 S) 
Inp —In(o) — 3° ap = Const; (15b) 


jedoch lasst sich T nicht mehr explizit darstellen, und die dem Poissonschen 
Gesetz, Inf — x Ing = const, entsprechende Beziehung (15b) hat ein gemisch- 


tes Zusatzglied P 
2 eeesSonn VE 
Pe <2 ae / ( a] 


erhalten. Die klassischen Adiabatengleichungen (7a) bis (7c) lassen sich indes 
formal weiterverwenden, wenn man — von den Differentialgleichungen (13) bis 


(15) ausgehend — 


oT oT 
peal CRESS Pe d(In Z,) a{T [o(InZ,)/OT lo} sha) 
dnt or atl oT lJ 


definiert und der Verlauf von x7 tiber T bzw. iiber plo fir ein betrachtetes 
_ Gasgemisch zahlenmassig bekannt ist. Fiir die Schallgeschwindigkeit im ru- 
henden Gas gilt dann beispielsweise weiterhin ‘streng: a2= %7 p/o, wahrend 
die Entspannungsbeziehungen (7a) bis (7c) mit Mittelwerten x7,7,'2u bentitzen 
sind, die im Iterationsverfahren bestimmt werden miissen. 

Die Stroémung mit veranderlichen spezifischen Warmen, aber noch vollig 
eingefrorenen chemischen Gleichgewichten, die von E. SANGER und dem Ver- 
fasser bereits in friiheren Arbeiten [4], [5] behandelt wurde, ist um so eher ver- 
_ wirklicht, je rascher eine Entspannung — von volligem Gleichgewicht im Kessel 
ausgehend — erfolgt. In der Raketentechnik ist sie bei Kurzdiisen am ehesten 
zu erwarten. Uberall da, wo diese Art von Stromung unter Nicht-Gleichge- 
_wichtsverhiltnissen auftritt, gilt beziiglich ihrer Entropie eine entsprechende 
_ Aussage, wie sie am Schluss von Abschnitt I fiir die véllig eingefrorene Stro- 


mung vermerkt war. 
Die Gleichungen (9) bis (12) liessen noch erkennen, dass die mittleren spezi- 


fischen Warmen C,or = H/T und C,or = U/T nicht mehr gleich den 6rtlichen 
_ Werten C,7 und Cyr sind. Man muss deshalb bei allen Stroémungsbeziehungen, 


_ insbesondere beispielsweise bei der Zeunerschen Gleichung fiir die Ausstrom- 
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geschwindigkeit wy aus einer Diise, den Unterschied zwischen mittleren und 
értlichen spezifischen Warmen sowie zwischen den aus ihnen abgeleiteten x- 
Werten beachten, der in den klassischen gasdynamischen Formeln verwischt | 
ist. Hierauf wurde bereits an anderer Stelle [6] ausfiihrlicher eingegangen. | 


III. Chemisch veranderliche isentrope Strémung 


Diese Art von Strémung, bei der sowohl die Veradnderlichkeit der Partial- 
druckverhaltnisse aller Gemischbestandteile als auch die Veranderlichkeit ihrer 
spezifischen Warmen beriicksichtigt werden miissen, tritt iiberall auf, wo Ener- | 
gieinhalt der inneren Freiheitsgrade und chemische Gleichgewichte sich unver- 
ziiglich den wechselnden Zustandsbedingungen anpassen kénnen. Sie erfordert 
die allgemeinste Form der Strémungsgleichungen, von der die in Abschnitt I 
und IT dargestellten Gesetzmiassigkeiten als Sonderfalle ableitbar sind. 

Wir beschranken uns in dieser grundsatzlichen Abhandlung der iibersicht-_ 
lichen Schreibweise zuliebe auf die Darstellung von homogenen Systeaey 
idealer Gase mit nur einer méglichen Reaktion, 


PAE Ay Was De Aa, 
=i t=1 


mit « Komponenten im energiedrmeren Zustand vor der Reaktion und w Kom- 
ponenten mach aufgenommener Reaktionswarme. Grundsatzlich lassen sich j je- 
doch alle weiteren Rechnungen auch durchfithren fiir Systeme von Gasen mit 
mehreren méglichen chemischen Reaktionen, seien es Molekildissoziationen, 
Rekombinationen oder Ionisationen, wenn man ausser den Molbriichen n,; der 
einzelnen Gase noch die Atomdruckbilanzen!) in die Rechnung einfiihrt. Die 
einzelne Reaktion sei gekennzeichnet durch: 

die Gleichgewichtskennzahl 


Ink, = S’y,; In b:;- 3», In Pr» (18a) 
i=1 bod 


die maximale Molzunahme (das heisst Molzunahme bei volligem Umsatz) 


Av= Sr, — Sm (8b) 
t=1 k=1 
die Reaktionsenergie?) 


AU=W,= = 2 U; — — 3% U,+ WW, (18c) 


1) Uber den Begriff «Atomdruck» siehe E, Scumrpr [7], S. 487. 


*) Hierin bedeutet Wo die nach erfolgter Reaktion bei 0°K yom Syste 
m pro v,, M;, aufge- 
nommene Reaktionswarme. i: i 2 ere 
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die Reaktionsenthalpie 

AH = W, =D)»; Hi — D3), Hy, + Wo = AU + Av ART, (18c) 

j=1 k=1 
die Reaktionsentropie 
o a Pe 4 ai 

AS = 2% Ae apa. S;, = ASo* + Ay oie = AUR ok (18d) 
die freie Energie der Reaktion 


AB = Slo, F,—))%F,= AU -TAS . (18e) 
i=1 k=1 


und die freie Enthalpie der Reaktion, 


AG = 3'1,6,—D/%, G, = AF + Av ART; (18f) 
4=1 k=1 


a 
3 worin >” jedesmal Summierung iiber die urspriinglichen Reaktionsteilnehmer, 
eS ) k=l 

» iiber die nach erfolgter vollstandiger Reaktion vorhandenen einzelnen Gase 
Bi=1 

bedeuten; y; kennzeichnet die Molzahlen der einzelnen Gaskomponenten nach 


- der Reaktion, »;, die Molzahlen der vor dem Umsatz vorhandenen Gase. Diese 
- Molzahlen stellen gleichzeitig die Héchstwerte (N;/Nz)may dar, die das betref- 
fende Einzelgas im reagierenden System erreichen kann. 

, Vor einer naiheren Untersuchung der Strémungsgleichungen unter den somit 
 beschriebenen Verhiltnissen ist noch zu bemerken, dass die Bedingungen kon- 
stanten Molvolumens dV = 0 und konstanter spezifischer Masse dg = 0 bei 

- verinderlicher Gemischzusammensetzung nicht mehr gleichbedeutend sind, 

' wie dies bei der Behandlung einheitlicher Gase oder eingefrorener Gemische 

~ der Fall war, weil nun Tetlchenzahl pro Volumeneinheit und Masse pro Volu- 

meneinheit nicht mehr einander konstant verhiltig sind. Es gilt vielmehr wegen 

; 


t V = f(m)/ge: dV_ _ dg {1 - £ } se} (19) 
0 (7a) 


fiir ein Gemisch, worin /(m,) eine spaiter genauer erorterte Funktion des Mol- 
bruchs 2, bezeichnet. Bei isentropischen Strémungsbetrachtungen interessiert 
die volumenspezifische Massenanderung, also die partiellen Ableitungen nach de 
- bzw. bei konstantem 9. In gleicher Weise versteht man unter Reaktionsprozes- 
~ sen bei konstantem Volumen solche, bei denen sich das Volumen pro Massen- 
- einheit nicht andert. 

| Betrachten wir zundchst wieder die allgemeine Gaszustandsgleichung (1). 
Sie hat bei chemischer Veranderlichkeit der Gase eine weitere Variable, das 


ve Q 
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z : . 

mittlere Molekulargewicht M= )3/n; M,, erhalten. Da dieser Ausdruck jedoch —: 

i=1 ; 

mit Hilfe des Massenwirkungsgesetzes und des van't-Hoffschen Gesetzes als: 
zusatzliche Bestimmungsgleichungen — dargestellt werden kann in Abhangig- | 
keit von nur je zwei der iibrigen, in der Zustandsgleichung vertretenen, veran- 
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Figur 1 
Linien konstanter mittlerer Molekulargewichte M (voll ausgezogen) und Linien konstanter Entro- 
pien S (punktiert) in Abhangigkeit von Druck und Temperatur, fiir Wasserstoffgas. 


derlichen Gréssen ~, T und g, so bleibt die Zahl der unabhdngigen Verander- 
lichen auf zwei beschrankt. 

Figur 1 zeigt die Veranderlichkeit von M in Abhangigkeit von # und T fiir 
ein reagierendes Gasgemisch vom Charakter X, = 2X am Zahlenbeispiel des 
Wasserstoffs. 


Die Veranderliche M kann wegen der Massenkonstanz 


t=1 k=1 


ee 
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Bich dargestellt werden in der Form M = My /(N/N_z), wobei Moy) = =p v, M;, 


das Gewicht der undissoziierten Ausgangsmasse, N/N; die ee Mol- 
zahl') der Gesamtmasse bedeuten. Die Molzahl N/N, lasst sich immer als Funk- 
tion F des Molbruchs 1, = N,/N eines geeignet ausgewahlten Einzelgases A des 
Gemisches darstellen, da bei einem System von z veranderlichen Komponenten 
die z Bestimmungsgleichungen fiir die z verschiedenen Komponenten, das heisst 
hier die z Molbriiche, immer auf eine Gleichung z-ten Grades fiir den Molbruch 
n, zuriickgefiihrt werden kénnen, dem alle tibrigen Molbriiche eindeutig zuge- 
ordnet sind. Allgemein gilt also: 


Wenn im gesamten Gassystem nur eine Reaktion méglich ist, lautet die 
Funktion F insbesondere: 


= F(n,) = —— (20a) 


“und man erhalt 


ee oe ee ee ee ee ee ee ee 
i) 


f(m,) = : , (20b) 
3 > 
4 k=1 
, 
wobei 7, alle Werte zwischen Null und val x y, annehmen kann. Es folgt weiter- 
hin: 
adf(m) _ _— tlm) (20) 
dn Vz as : 
; Avy 4 


| Mittels dieser Beziehungen kann man nun die allgemeine Gaszustandsglei- 


chung in der Form 
i p= eke tf (20) 


schreiben. Ihre partiellen Ableitungen lauten dann: 


T[0n,/0T}y 


s<l, =+ . {r Ne 1} + oa tp a i}, (21) 


eee Spam 


1) N bedeutet die Zahl aller jeweils vorhandenen Gasteilchen, Na die Zahl der jeweils vorhan- 
- denen Teilchen der Gasart A, bei konstant bleibender Gesamtmasse. 
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ee i eee 
lm Ele Nee a= 8 at 
aerlpsctigelPa restanea or mney As ae 
und 
alata ee Un aioe an 
Sie gehen fiir 
dj(n,) = ee] te, [ ee] La [ “| = aed [ aa |= [ ae 
p Te 2 
in die klassischen Grenzfalle 
Ea ey a i 
a] rT 19 a ; 
itl +£ 
al rT us - : 
a4 ae z 
und 
Eanes 


ber. Aus den Gleichungen (21) bis (26) folgt ferner: 


ap _ ag 
P ) 


af (m4) 
f(m,) ° 


|-° 


(21a) 
(22a) 
(23) 
(24) 


(25a) 


(26a) 


(27) 


Fiir Reaktionen mit unveranderter Molzahl Ay = 0 kann die Zustandsgleichung 
wie bei einem chemisch unverdnderlichen Gasgemisch behandelt werden. 

Ahnlich wie die Gaszustandsgleichung lasst sich auch die Entropiegleichung 
(4) als Funktion der vier Veranderlichen p, T, 0 und m, auffassen: 


s = Di piln,) yi 
1=1 


(28) 
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worin wieder 


p iy Pi 
S.= soe4 [c, > — Arf? (28a) 
0 0 
Q a Ad 
= see4 fc, 3 - AR fe (28b) 
0 0 


mit molaren spezifischen Warmen entsprechend den Gleichungen (11) und (12), 
und worin ferner 


cay te Di (28 ) 
V:(M) = Om os = 7 A c 
A fee 
Mo (1 v, n,) 


Wenn A eine der m Gasarten darstellt, die sich nach erfolgter Reaktion bilden, 
so gilt fiir alle m; zwischen (z — w) und z: 


- wahrend fiir die Molbriiche der restlichen « Gasarten, die wir zwecks Unter- 
 scheidung mit der Laufzahl k versehen haben, 


Vy i=1 
Ny = — 1 moa 
A 
Dk 
k=1 


- geschrieben werden muss. 
Nun lauten die Ableitungen der Entropiegleichung, in Funktionen von 


(p, T), (e, T) oder (f, @) ausgedriickt : 


BP tan ZAP 9 P54 ooo [set] 7 | 
= (29) 

£ (Re), 5+ oes [285], ] 48 | 

ds, p= y (| ee) Sit Palm) aad a | 
a (30) 

a S[eee], Set oat [Fe] 42: | 


ZAMP VI1/4 


} 


50) IRENE SANGER-BREDT ZAMP- 
bzw | 
is Ss (ree ] Sit Vil) =I, dg | i 

(31) 

ye ee 


Hierin sind die partiellen Ableitungen der molaren Entropien einzelner 
Gemischkomponenten, [0S;/0x],, weiterhin gemass den Ausfiihrungen in Ab- 
schnitt II fiir reine Gase mit verdnderlichen spezifischen Warmen zu behandeln. 
Man erhalt entsprechend: | 


[s=I, LAR {| 0(InZ,) |,+ O{ T[O(InZ,) /0T}o} |,+ 1} | 


oT oT oT 
a (29a) 
a | a 
ae Ane ee -=8; (29) 
=I, oe: {| — ],+ [ ete oleae 1, 
a (30a) 
che 
boleh ast ibe ers (30b) 
Galan (31a) 
und 
3], BH, - +S av 


Jeder der Ausdriicke [0q; (m,)/0x], in den Gleichungen (29) bis (31) lasst 
sich auflésen in ein Produkt: 


Apna) 2] 
y? 


dn, dx |; 
wobeli 
ae Bags 
aQ,(n,) _ Pi (M4) "A RoI 
ati i (ore Ag. \a (32a) 
a Ny Moo (1 ~<"n, ) 


othe ® the 


ye ee ey ee ee a ee ee ee a 
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bzw. 


Vygtee 
al Vy 
ap y,(%4) CA) ~ : (32) 
an, Mog (1 ‘a aay 
2 


wird. Von den moglichen Varianten des zweiten Faktors [07,/0x], in obigem 
_ Produkt kann man die isothermen Anderungen des Molbruchs 1, mit Hilfe des 
_ Massenwirkungsgesetzes (18a) — in der Schreibweise: 


E,(n,) = K,p-” (33a) 
oder 
&,(1,) = Fis Crear, ee ? (33b) 
worin 


yp 
G1 R=1 ot 
k=l 
= €, (n,) [F@))”, =. 
erfassen. Dann wird: 
Op Ir [0€,(m)/dnj]r =» Eg (m,)/dn, pp 
und 
4] 2 [0& (73) /Oolr Eo(7) Ay : (35) 
Og Ir [0E,(m)/dn,]r dE,(n,)/dn, © 


Die Molbruchanderung in Abhangigkeit von der Temperatur wird unter zu- 


_ sdtzlicher Anwendung des van’t-Hoffschen Gesetzes: 


OK, > AH 
ie |= P ART? 
oder 
OK, AU 
Ga ],= Ke ART?’ 


_ jenachdem Druck oder Dichte wahrend der Anderung konstant bleiben sollen, — 
_ bestimmt zu: 


oy aE aaa +p am ART (36) 
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und 


[7], = Let], [or]. = + eee RTF ie 


In den Gleichungen (34) bis (37) betragt dabei 


at SO EMG BMS 9c nee Se (38a) 
ae oe V7, PP a | 
4=1 
bzw. 

o bye 

dé, V; Et Ay 
oe = &,(n,) Pits: Es o zs v,/Av— nN; (38b) 

A t=1 ea| 3 — Ny | 
i=1 


Die beiden restlichen partiellen Ableitungen von m, endlich lassen sich mit 
Hilfe der Gleichungen (25) und (26) auf die Beziehungen: 


243 ey, Ae. (39) 
bzw. 

On) _ 1 a a 

Le |, = i 1 7] [earns ie 


zuriickfithren, in denen [0n,/0T], und [0n,/0T], aus den Gleichungen (36) und 
(37) bekannt sind. 

Man erhalt so — unter Beniitzung der Gleichungen (38) bis (40) — fiir isen- 
tropische Entspannung eines chemisch verdnderlichen, homogenen Systems 
idealer Gase mit nur einer méglichen Reaktion und mit veranderlichen spezi- 
fischen Warmen die folgenden Differentialgleichungen. 

1. Fiir die Beziehungen zwischen Druck und Temperatur gilt: 


(dg; (: 0 OS; 
ar Onn [ort], Se+ eed) [Gr'],} a7 
(4 


+3 ee [Set] Set pale) self a = 


t=1 
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i 1S AH AS k=1 i=1 dT 
ART Mo (va — Av n,) T 


(03) a o4 gz 
| (»: ee n,) v4 Ny pee Vie VA VED, 1, Cy, 
ce ~ 


ds, p= ~ 
Mo (¥, — Av mi)? dy % 
41= 


@o og [e4 
(m— Srms) am dy Sm AR 
i=1 k=1 k=1 dp 
Ay AS + Vi 


| W 
| Mo (v4 — Av ny)? Dd); oo (Ya yn) | Pp 
i=l 


wo lod 
(».- D5 m,) Van, Av > % 
i 1 


= AH (AH — AG) 


ART? 


t=1 


(3) 
Mo (va — Av mg)® 3% 
=1 


io z 


4d) Ye Dy Mi Cp, 
k i=1 


=1 
= Moo (va — Av 14) if 


t= 


@o a 
= ; A 
(>, 2 V; n,) V2 14 Pe: Av (AH — AG) 


eet 
= 
— 
2 
| 
& 
= 
= 
Yo 
Me 
Ag 
isl 


Hierin betriagt die Enthalpie pro Kilogramm: 
(> ats m) v2 24 Av 2% AH 


2 i dp (41a) 
Moyo {vz — Av na)? Dd)»; 
t=1 


dh = c, aT 


mit 
a 4 w \ (o4 
> Me Dy Mi Cn, (4-3 V; m) VAN, ane Vip ote 
k=1 i=1 tee ele Sein 
o ART? 


Moo (va— AY n3)> D257 
t=1 


ay M qo (¥, — Av mj) 


1) Wegen H—G=TS. 


| 
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und ferner die reversible Arbeit pro Kilogramm: 


dA 
M M9 (v, — Av m4) Z 


tei ae dp 


(>, hier ta n,) vm, Av S) vy, | 
4 j=1 ey fe ( AH aT a #) 


wo 


M oq (v4 — Ay n,)> d'v, 


3 | 


Da der Anteil der Reaktionsarbeit an der reversiblen Gesamtarbeit unter 
anderm auch temperaturabhingig ist, gilt bei chemisch verdnderlichen Syste- 


men nicht mehr: 
OH OS 
Gok ae Dae 


Nachdem man iiberall durch das mittlere Molekulargewicht 


ot 
4 DM 
k=1 


Mas 
M 9 (¥% — Av n,) 


gekiirzt hat, lasst sich Gleichung (41) auch in der Form 


(.- oo x» ; na Ap 


AH AS : 
a ‘ARE ane Cy, 
(vy— 1, Av)? Y'y, <é: 
Cee a ‘ dT 
: : iT (414) 
(u—m 3») mn, Ap 
ee AvAS+AR 


wo 


(vy — _, Av)? Sy, 


t=1 


schreiben, und ist so unmittelbar mit den Formeln (6a) in Kapitel I bzw. 
(13a) in Kapitel II vergleichbar. 


2. Fir die Beziehungen zwischen Dichte und Temperatur gilt: 


Oe 2d {2a [sz], Sit y,(n,) sz. ar | 
{ dp,(n,) Ton 0s a 
SLI [Fe], 5 eg [954 a= | 


} 


Cee eR A RN Pee a ee ee ee 
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wm ied 
(1m 3» vm, dv 3 
t=1 =1 


ds, 7p = wo (a) 
| (vy — ny Ayr) J(os— man? Mt (v= m1 Sw.) ny APy| 
t=1 t=1 
og 2 
v v n,C 
AU AS 124 ee : «| aT 
ART Moo (v, — 2, Ar) | ah 


t=1 


M oo (v4 — 24 AP) (ox A)? Sand (v,—m, 3 ») ny A] 


(4 
, >) v,AR 
Ae a 
M qo (va — 23 A”) | & 


t=1 


o e4 
(m= ma, ») m1 Aw DM 


Mo (¥, — 1% AY) [tos — ma A? Du + (v.—m x» n, A>] 
4=1 t=1 
a & 
v y N;,Cy. 
AU(AU— AF) | a aire oem «| dT 
Ase D2 Moo (v4 — % A?) | FE 


@ 
| (4m 3%) vim Ay 2% 
i=1 


te=1 


Moo (va —- 1 AY) J Nn Av)® S)% “+ (mu—m, Ss) md] 


a 
Vz HAR 

Av(AU — AF) k=1 
a Mo (v%,— 2 AY) [ @ 


55 
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Bei der Ableitung dieses Ausdruckes wurde die Beziehung U — F = TS be- 
ntitzt. In Gleichung (42) betragt: 


ow a 
s y 
(», hy “) Vv, Nn} Av Vy, . 
q=1 k=1 . 


du=c, a1 = = 
M oo (¥, — 2, Av) [Co A)? + 14+ (a—m x») nA] | 
- t=1 i=1 ( 
“i 
x —\d 
0 | e 
mit 
a & 
VA Dy, Den; Cy, 
“* pl fa 
*o Moo (74 — 4 40) 
wo (od : (4 
5 P A ; : 
si (2) a) A hT=> 
M yo (v2, — 2 Av) J(r— 1, 40)? Sine + (=m Ss») nz a%| 
= = 
Ferner ist: 
a 
v v, AR 
Be ‘a aoe k do 
M as cu @ 
@ a 
(.—m, S| vn, Avy Dv, 
mal k=1 
| = (42c) 


wm oO 
M oo (¥4 — 2g Ar) — 1, Av)? Sy, + (>, — nN; x») Nz ary| 
t=1 i=1 


Im iibrigen gilt in Analogie zu dem bereits unter 1 Erwahnten, dass bei chemisch 
veranderlichen Systemen wegen der Temperaturabhangigkeit der reversiblen 


Reaktionsarbeit : 
( as aes 
OT }» Ge i 


wird. 


ee ee Ce ee ee a ee a 
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In der Form: 


@ 
(v.— ms 3m n, Av 
i= 


ANOS ASS - 
@ @ SS ARP en 1, Cy, 
(v4 — 14 Av)? D'y, + (vm, x») n, A*y a 
_— i=1 t=1 dT 
a : : (42d) 
(».- nm 3%) Nz Ay 
as Ay AS+AR 


(vy, — 1, Av)? D'y, + (v1— 1, 5.) n, A*» 
ASI i=1 


ist die Beziehung (42) wieder ohne weiteres mit den entsprechenden Ausdriicken 
(6b) in Kapitel I bzw. (14a) in Kapitel II vergleichbar. 
3. Fiir die Beziehungen zwischen Druck und Dichte bei isentropischer Zu- 


standsanderung ergibt sich aus den Gleichungen (41d) und (42d) durch Elimi- 
nation von dT /T der Ausdruck: 
£1 de;(n;) [222] as; 
£ S;+ 9,(n [ ‘| 
dp ie 2 dn, 00 p P;( a) 00 p @ do | 
gr 5 { 2p; (%2) [] Sheets 
a S; } (n i] 
24 dn, ap ig¥:) = nile Op Io) 7 
v,— nN v;,) 1, A 
(» ot : ARATE SS ole 
2 Vane te 
(vz — mq Av)? Dd*y,; 
a i=1 
(» — ity ») Ny Av 
i=l AU AS id 
as a arr + Am Su | (43) 
(v, — nm, Av)? n+ (4 —m 3) n, A*y 
i=l i=1 
(=m, yn) Ny, A*y AS 
: 4=1 Tf AR 
(v, — ny Av)*® D'y;, + (4 m 3%) n, A*y 
w t=1 B t=1 2 ’ dg 
(n—m 3) n, A*y AS e 
’ 4=1 = a te AR 
(vy — 1, Av)? D'»; | 


t=1 


~ mit dessen Hilfe man gemiss der Beziehung dp/do = a” unmittelbar die Schall- 


geschwindigkeit im ruhenden, chemisch veradnderlichen Gasgemisch bestimmen 
kann. Gleichung (43) entspricht den Ausdriicken (6c) in Kapitel I bzw. (15a) 


in Kapitel II. 
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Die Gleichungen (41d), (42d) und (43) sind nicht mehr ohne weiteres inte- 
grierbar. Wenn man bei endlich grossen Zustandsanderungen in chemisch ver- 
anderlichen Gasgemischen die klassischen Adiabatengleichungen (7a) bis (7c) 
formal weiter anwenden will, so hat dies mit Mittelwerten zu geschehen, die im 
Iterationsverfahren bestimmt werden kénnen, wenn der Verlauf der 6rtlichen 
Werte im Koordinatensystem der Zustandsgréssen , 9 und T zahlenmassig 
bekannt ist. 


Ausserdem sieht man es den Gleichungen (41d), (42d) und (43) ohne wei- | 
teres an, dass Gase, die bei unveranderter Molzahl reagieren (das heisst Ay=0), © 


hinsichtlich aller Entropiebetrachtungen wie chemisch inerte Gasgemische im 


Sinne des Abschnitts II behandelt werden kénnen, das heisst wie Gasgemische, 


fiir die dn, > 0. 


Definiert man die Grésse ‘ e@ ap 
POLE: 


als den differentiellen Isentropenexponenten — analog zur Grésse x im Poisson- — 


schen Gesetz — und halt man die Analogie zu den klassischen Adiabaten- 
gleichungen auch fiir die Beziehungen zwischen Druck und Temperatur sowie 
zwischen Dichte und Temperatur aufrecht, indem man 

ti pen an Ope: 

Pp Cw Op, G eas i 


und 
Ladones of 
oar ~ am —1 


setzt, so ergibt sich: 


v,—N "») Ay | 
ke SHAS ie 
@ ZET ts 
(vy, — , Av) Dy, e 
23s t=1 
0 23) 
V4 — Ny, *) n, Ay 
BS ( a AU AS : 
: o 0) ART aot a ee 
(va — 1, Av)? S'y, + (m1 —m Sy n, A*y oa 
t=1 t=1 


2 = -+ AR 
(v,— Av)? Soy, + (s4—m | Ny A*y 
~ i=1 i=1 


(v1—m, 3) n, A®*y AS 

=______ + AR 
(vz — N14 Ay)? PIE & 
t=1 


- t (44) 
(1. Ny Sy) n, A*y AS 
t=1 


ay 


eo. leit) Rh bh ee PORE CeO NT ee ee 
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i 0 ART toon 
srl T (vg — ny Av)? S; = 
Co a oT P t=1 
pT Os Os o , (45) 
ie p 
qoOLE p l Op T n1— m 3m) n, Ay 
=A AU AS a 
2) PAR onde + Cy 
(v, — m, Av)? Dy, ech 
1—1 
bzw 
Os Os 
tah Pane 
OoT re Os 
soln 
( o 
v,— nN v;) m,Av 
aay : : AAS 
: b ART + os % Op (46) 
(vz — ng Av)? Dv, + (vs — m3) nN, A*y ag: 
t—1 v=1 
v,~—n v,) nm, Av 
(» om ‘ 3 AAS 4 ee 
@ @ a Bea + Dn; Cy, 
(v, — 23 Av)? DS’y, + (= , 3%) n, A®y te | 
tad tot 
Es ist 


o> 1 
a a Op,T oO me 


i A 


_ Die Isentropenexponenten 04,7, 0),7 und oy, fiir die Beziehungen zwischen je 
- zwei Zustandsgréssen sind also nicht untereinander gleich. Insbesondere sind 
- sie auch verschieden von dem Verhiltnis der 6rtlichen Gemischwaérmen 


[7 
a OT |p 
am |) 
OL ip 
wm 
(.-m 3») n, Av an F 
1=1 
a ARTs tem Cn, (47) 
(v,— mg Av)? Dd, »; 
41=1 
—n n, Av 
(» no) ») : AeU : 


w 


(vy, — 2, Av)? is (>, — m3») ny Ay 
4=1 


t=1 
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Fiir dn, -> 0 streben jedoch alle demselben Grenzwert 
De nN; Cy, 
i=1 


Zz 
t=1 


Zu. 
Unter Verwendung der bisher gewonnenen Erkenntnisse erhalt man fiir das | 
spezielle Beispiel eines nach dem Schema X,= 2 X reagierenden Gases, wie 
Wasserstoff, also mit 


tO (DE ke Di vga 2s eipe= 1, Av=1,  %1 = ty =2Z; 
5 g=1 


m,=Ny, My =My,, AH=AU+ART=2H,+ Wy — Hy,, 
AS =2 Sy—Sy,: 


Ny(1—ny) AHAS 


o,,= (2— my)? ART b Mr Coy +  — te) Coy 
P,0 iS AC AS 
Jong ARDS Roe ee a oe, 
(44a) 
et 3 ty ‘ 
Ny (1 — ny) 
Go apen 4 tae 
Ny(l—nyz) AHAS 
O,7p= (2—ny)* ART tether Gage ed ta Con, 
pT ny(1— my) AUAS ; (49a) 


—— ART + Ny Co =>- (1 = Nz) Con 


ny (1— my) AH AS 

C7 = (4 — 3 ny) “Ger ot tH Coy + (1 — My) Con. 
(i ny (1 — ny) AU AS sf (46a) 

(4—ingi~ ART * Cog? Me Cup 


und 


Mz (1 — 177) AlMerk 
ee Te le “SA Rep ee eee ea 
ts a ae oo ae 
(4— Shy) ARTE + MH Cogt (1 MH) Cog, 


é In den Figuren 2 bis 5 sind in Abhangigkeit von Druck und Temperatur 
die drei verschiedenen Isentropenexponenten %p,r, %,r und o,, sowie das 
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Figur 2 
Linien konstanter drtlicher Isentropenexponenten & ’,T (voll ausgezogen), Linien konstanter ort- 


licher Isentropenexponenten ©, 7 (soweit sie vono »,T abweichen: gestrichelt) und Linien konstanter 
Entropien S (punktiert) in Abhangigkeit von Druck und Temperatur, fiir Wasserstoffgas. 


Verhiltnis der értlichen spezifischen Gemischwarmen xr und das Verhaltnis 
der mittleren spezifischen Gemischwarmen 


& 


» % A/T 


t=1 


Dm UilT 


t=1 


Co 


fiir dissoziierendes Wasserstoffgas eingetragen, um die Unterschiede zahlen- 
massig zu veranschaulichen. Punktiert wurden ausserdem in allen Abbildungen 
Linien konstanter Entropien gezeichnet, langs denen Zustandsanderungen ent- 
sprechend den Voraussetzungen dieses Abschnitts zu erfolgen hatten. 

Die dargestellten 6rtlichen Isentropenexponenten zeigen zunachst die zum 
Teil erheblichen Abweichungen untereinander sowie vom Verhaltnis der spe- 


62 TRENE SANGER-BREDT ZAMP 


3 NY F gRSAgsswegssaggss 2 S S 
5500" SNC 
5000 178 
4600 Hit +s TET ios a Si : i 
Hof ‘ 2 ll 
4000 ak 
S A “f T 
pasa0 18 =. cml 
ce “HI 52 
x eh WA | L| TS 
q 122 5 - - = 
83000 : an = APL pee 
& ii S HLL | ZO TAI i 
N 2500 128 Bf TPs 
130 sy. 
2000 V2 
hh tt 
soo NTT Tt ici 
1000 498 HH 
Pes ul =a Paes! : ee é 
500 fH NT Tc TT THbo by 
seed ere = ae nar . Se el ae 0 
id ain seen 6 
700 10 7 Gi Druck plat] Goi GU 
Figur 3 


Linien konstanter 6rtlicher Isentropenexponenten Oy o (voll ausgezogen) und Linien konstanter 
Entropien S (punktiert) in Abhangigkeit von Druck und Temperatur, fiir Wasserstoffgas. 


zifischen Warmen bei konstantem Druck und bei konstantem Volumen im Fall 
eines chemisch verdnderlichen Gasgemisches mit Gleichgewichtseinstellung. 
Wegen ihrer starken Veradnderlichkeit iiber p und T sind sie von praktischem 
Wert vor allem bei der Berechnung von Zustandsainderungen innerhalb enger 
Bereiche von Druck und Temperatur unter Verwendung der tiblithen Stré- . 
mungsbeziehungen — besonders also zur Bestimmung der Schallgeschwindig- 
keit im ruhenden Gas bei tragheitsloser Gleichgewichtseinstellung beispiels- 
weise. Auch bei zwar endlich grossen, aber noch geringen Zustandsanderungen 
kann man mit Hilfe derartiger Auftragungen, wie sie hier fiir Wasserstoffgas 
vorliegen, fiir irgendein chemisches Gleichgewichtssystem die klassischen, sehr 
ubersichtlichen Strémungsbeziehungen weiter anwenden, wenn man fiir x je- 
weils die richtige Beziehung einsetzt und aus dem betreffenden Diagramm 
dafiir einen Mittelwert langs der betrachteten Zustandsanderung bildet, ahnlich 
wie dies bereits in Abschnitt II vorgeschlagen wurde. Bei Zustandsainderungen 
tiber weitere Druck- und Temperaturbereiche kann man mit geeignet gewahlten 


_— >) 2 
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Linien konstanter 6rtlicher Adiabatenexponenten x (voll ausgezogen) und Linien konstanter 
Entropien S (punktiert) in Abhangigkeit von Druck und Temperatur, fiir Wasserstoffgas. 


Mittelwerten zwar noch zu Ndaherungsergebnissen gelangen; doch ist bei sehr 
weiten Anderungsbereichen meist die Anwendung von kombiniert graphisch- 


_rechnerischen Verfahren unter Verwendung des Entropie-Enthalpie-Diagramms 


dem Versuch vorzuziehen, auf die hier dargestellte Art auch Tafeln fiir muttlere 
Isentropenexponenten zu entwerfen, die immer nur sehr spezielle Verwendungs- 
fahigkeit haben kénnten, da sie nicht nur stoffabhangig sind und von der Be- 


_zugstemperatur der dargestellten Temperaturspannen abhangen, sondern auch 


von den Bezugsdrucken der dargestellten Druckdifferenzen. 


IV. Anwendungsbeispiele 


Zam Schluss soll noch am Beispiel einiger besonders haufig benutzter Stré- 
mungsbeziehungen gezeigt werden, welche Werte an Stelle des klassischen x 
treten miissten, wenn die Beziehungen auch in den Fallen chemisch verander- 
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Figur 5 
Linien konstanter mittlerer Adiabatenexponenten %o 7 (voll ausgezogen) und Linien konstanter 
Entropien S (punktiert) in Abhangigkeit von Druck und Temperatur, fiir Wasserstoffgas. 


licher Gleichgewichtsstrémung anwendbar bleiben sollen. In der klassischen 
Gleichung fiir die Schallgeschwindigkeit zum Beispiel miisste, wie bereits er- 
wahnt wurde, x durch den Grtlichen Wert o,, ersetzt werden: a? = Op,o Plo. 


Die Zeunersche Gleichung fiir die Ausstromgeschwindigkeit W » aus einer 
Diise mit linearen Strémungsverhiltnissen: 


=> ieee Ee Go lone 


ware demnach beispielsweise in der Form: 


x F = 1 & Z 
Wm = \/2 = OT, e Po 1 “OT m b ly ($2) (6 m0 m GAC 


“or, “ot, TAN Pm 
(48a) 


zu schreiben, worin der Index 0 die Zustinde im Ofen, der Index m die in der 
Diisenmiindung bezeichnen. 


Oe 


vTseqe 


el el 
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Die Bestimmungsgleichung fiir die Geschwindigkeit w, im Diisenhals einer 
Lavaldiise: . 


an fob eR RO “ 


v 


wurde 


(49a) 


“OT y é Po ae ee ete 
oT, ah Qo Py» 


bei Beriicksichtigung der verschiedenen x-Werte lauten miissen. Hier bezeich- 
net der Index v iiberall die Zustande im engsten Diisenquerschnitt. 
Die durch Gleichsetzen der Geschwindigkeiten in (48) und (49) sich erge- 


bende klassische Beziehung fiir das Verhaltnis der Drucke im Ofen und im 


engsten Querschnitt: 


Pe $4 C : Avene ‘ (50) 


miisste bei Ableitung aus den Gleichungen (48a) und (49a) genauer 


(50a) 


Pore (tose ia meter eA OT 1 jp dale 
—1 


2 Mae OT 
heissen. 

Die klassische Gleichung endlich fiir das Verhaltnis von engstem Quer- 
schnitt f, zu Diisenmiindungsquerschnitt /,,,: 


fe e cae io eed ee E “ alate (51) 


miuisste genauer 


te [ee st tom | Fb, 0% Fl x Tse * Le 1 Maen eo eom)! (Po eon 2) 
ita i Po St 

Fi - Le ENTE Biehl pte) 
x |/ Tp, 04 E 0Tm “OTs ay Pom: 20m ) pork) 


= x —i1 
“by, 00 . sd Olen Po 
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geschrieben werden. 
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Summary 


The notion of an adiabatic exponent—as the ratio x of the specific heats at 
constant pressure and at constant volume—which had been so productive in clas- 


‘E. SANGER und I. Brept, Uber einen Raketenantrieb fiir Fernbomber, Dtsche. 
Luftfahrtforsch. UM 3538, Ainring 1944 (Edition Robert Cornog, St. Barbara, | 


; 


sical gasdynamics, gets unfit for thermodynamic flow, if the suppositions of con-_ 
stant specific heats and of constant partial pressures can no longer be maintained _ 


for the different gas components. 


In these cases, the heterogeneous origin of the different, only apparently uni- — 


form exponents, which occur in the flow equations, becomes evident and their 
numeric values begin to differ from each other. For dissociating gases, e.g., the 
adiabatic exponent is no longer identical with the ratio of specific heats. 


In the present paper, the influences of variable specific heats and partial : 


pressures on the adiabatic exponent are taken separately into account. With new, 
appropriately determined exponents o, generalized equations for isentropic flow 
are established, from which the classical formula can be derived as a limit case. 


(Eingegangen: 8. Februar 1954.) 
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Zum Stromrauschen von Halbleitern 


Von WILLY BAUMGARTNER und Hans ULRICH THOMA, Ziirich?) 


Verschiedene wichtige Ergebnisse der Rauschmessungen an Halbleitern haben 
bis jetzt noch keine Erklarung gefunden, die auf Ad-hoc-Mechanismen verzichten 
kann. So ist die Grosse des Rauschquadrates nicht verstandlich, ebenfalls nicht 
seine Frequenzabhangigkeit (zum Beispiel f~1-Spektrum, das bis zu tiefen Fre- 
quenzen lauft)?). Nimmt man an, dass speziell fiir letztere das Spektrum unkor- 
reliert aufeinanderfolgender Einzelereignisse massgebend ist, so muss schon das 
Einzelereignis selbst eine grosse Ablaufdauer aufweisen. Die sonst plausible Gi- 


solfsche Theorie*) scheitert in dieser Hinsicht an der zu kurzen Dauer der freien 
Beweglichkeit der Elektronen und Lécher. 


1) Abteilung fiir industrielle Forschung des Instituts fiir technische Physik der ETH. 

2) A. VAN DER ZIEL, Physica 16, 859 (1950). K. W. Borr, Ann. Phys., 6. Folge, 14, 87 (1954), 
D. K.C. McDonatp, Rep. Prog. Phys. 12, 56 (1949) 

3) J. H. Gisorr, Physica 15, 825 (1949) 
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Es liegt nun nahe, den in der Theorie der Photohalbleiter und Phosphore wohl- 
bekannten Vorgang des Einfangens in Haftstellen zu beriicksichtigen. Das Einzel- 
ereignis hatte dann folgenden Verlauf: 

1. Bildung eines Elektron-Loch-Paares (thermisch oder optisch). 

2. Wandern des freien Elektrons (das Loch habe vorerst die Beweglichkeit 0) 

wahrend einer gewissen Zeitdauer. 

3. Das Elektron wird eingefangen. In diesem Zustand sll keine Rekombination 
_ mit Loéchern stattfinden. 

4. Reemission aus dem eingefangenen Zustand und Weiterwandern des Elektrons. 

5. Vorgang 3 und 4 mehrfach wiederholt. 

6. Rekombination des freien Elektrons mit einem Loch. 

Da die Verweilzeit im eingefangenen Zustand bekanntlich grosse Werte er- 
reichen kann (= 1s) so besitzt das oben beschriebene Einzelereignis die geforderte 


_ langfristige Dauer. 


Die Rechnung ergibt fiir den besonders interessierenden niederfrequenten 
Anteil des Stromrauschens: 


Ww ADT pS er gt 1 : 
L* Bre 1+ 0% a? (1 + @/f)? 
Dabei ist 
W = Powerspektrum, das heisst - W df = 1?; 
- N,= Anzahl freier Elektronen i Parallelepiped mit Lange L cm; 
wo = Kreisfrequenz; 
e = Elementarladung; 
E = Feldstarke; 
, = Beweglichkeit des Elektrons; 
a = mittlere Lebensdauer des Elektrons in den Haftstellen «> £, 0; 
B = mittlere Lebensdauer des Elektrons im Leitungsband, falls keine Rekom- 
bination auftritt; 
0 = mittlere Lebensdauer des Elektrons im Leitungsband, falls kein Einfangen 


auftritt. 
Fiir die entsprechende Rauschziffer FR ergibt sich 


Bee ae 10m ja ba tite 
R= kT B+e 1+ 070? (+ 0/6) 
mit 
ee L0H, op = 10-*s,. fi 10-45; 
rF =40V-1, pw, = 10? mn ise. == LOA Wem: 
folgt 
i R,-9~ 40- 


Wir erhalten also wohl die bis zu tiefen Frequenzen gehende Anderung des Rau- 
schens, aber nur massig grosse Werte der Rauschziffer. Haben aber die Locher 
eine Beweglichkeit p,, +0, ist also wahrend der Haftdauer eines Elektrons 
laufend ein Stromtransport durch entsprechende (Ladungsneutralitat!) Locher 


méglich, so folgt als Rauschziffer: 


ee eine ye. 1 = 
aS RE B  1+o2a*(1+e/6) 
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was mit den oben angefiihrten Daten ergibt (mu, = 102 cm?/Vs) 
eg sare 


(2) 
also sehr grosse Werte. 


Selbstverstandlich kann eine Ladungsneutralisierung auch durch ein in der — 
Nahe des Elektrons ebenfalls festsitzendes Loch erfolgen, so dass dieser Wert — 


nur als obere Schranke zu lesen ist. f 

Auf die einleuchtenden Beziehungen mit Ankling- und Abklingzeiten des 
Photostromes in Photohalbleitern!) und der Phosphoreszenz von Luminophoren 
sei hingewiesen. So besteht die Méglichkeit, durch Annahme einer kontinuierlichen 
Haftstellenverteilung, wie sie GARLICK?) zur Deutung des Abklingens von ZnS-Cu- 
Phosphoren verwendet, eine Frequenzabhangigkeit proportional f-* (1 < k < 2) 
zu erklaren. 


Summary 


The wellknown trapping of electrons and holes seems to give an explanation 
of the current noise in semiconductors. It is pointed out how especially the 
magnitude and low frequency dependence of the noise power can be understood; 
the possibility of a strong correlation between low frequency noise and photo- 
conductive properties of semiconductors is emphasized. 


(Eingegangen: 10. Oktober 1954.) 


Vibration of Isosceles Triangular Plates 


By Hueu L. Cox’) and Bertram KLern$), Los Angeles, California 


1. Introduction 


Frequently in the analysis of structures it is necessary to determine the 
lowest natural frequency of flexural vibration of a triangular plate. The approxi- 
mate solution presented in this paper is based essentially upon the method of 
collocation, i.e. for an assumed deflection function, which satisfies the boundary 
conditions, the governing differential equation for the plate is satisfied at a finite 
number of points. Comparisons are made for certain cases with other known 
solutions, and the agreement appears to be satisfactory. 


2. Simply Supported Plates 


Let # and y be coordinates in the middle surface of the plate as shown in the 
lower part of Figure 1. The static differential equation for a uniform plate sub- 
jected to a uniformly distributed lateral load is 


DY" yp = 7, 


1) Uber einen experimentellen Befund in diesem Sinne vgl. R. L. Perrirz, Proc. IRE 40, 1440 
(1952), speziell S. 1451. 
7 


G. F. J. Garrick, Luminescent Materials (Clarendon Press, Oxford 1949), S. 33. 
8) Research Engineer, T, Y. Lin and Associates, 


ae 
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Figure 1 
Vibration Coefficients for Triangular Plates. 


plate stiffness ; 


w = deflection of plate, positive downward; 
q = lateral load per unit length squared. 


69 


If the term g in the above equation is replaced by an equivalent inertia load when 
the plate is vibrating, there results 


where 


2 
DPEP3w Aw, 


@ = natural frequency; 


gravitational acceleration. 


The boundary conditions are 


( 


02w 


! 


(W)y 4 = ()z-4(almyy = 9» 


Oy? 


i 


0?w rt (<> 4° ar) 
“! Ot? ae On® Or? x=-+ (a/h)y 


(1) 
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where 7 is the normal direction to the lines 


a 
Ope res 


and ¢ is the tangential direction of any line. Along a rectilinear edge the term 


0?w/dt? is equal to zero. 
The deflection function 


a Pi 
w= (4,7, + «f+ «3 fF) cosk ear: we (4) 
where 
Sy oe) ee ee ae Bh OAS ian tag ine Soy fo) eee, 
A= 2sin* 75 sin; By 2 Si 2h 82? Fa —=-2/ Sin 7 a a 


« = generalized coefficient; k = odd integer 

satisfies equation (2) exactly and the condition 
( 0?w = 
On® Jx=+ (a/h) y _ 


from equation (3) at some point in the region h/2 < y < 3 h/4; the condition that 


0?w 
=> = 0 
| Oy" es 


is satisfied exactly only when + = 0. Selected points for collocation will be on 
the line ¥ = 0; therefore, no serious error should result from not satisfying 
completely the boundary conditions. The first term of the deflection function 
given by equation (4) appears physically to have the shape that one might expect 
for a triangular plate vibrating in its fundamental mode. 

It will be necessary to differentiate the deflection function given in equation 
(4) and substitute the proper derivatives into equation (1). 


Let 
Fiyj=y4 Fi +%e%R+ oF. (5) 
Equation (4) may be written as 
Ree PE 
phys Oe ee ee (6) 
Derivatives of w at x = 0 are: 
(%),-9 = F(y) , 
02h XR h 2 
Cree aia Re ” 
04w rt h- \4 
(it),o7 FO [QS EY, (8) 
0°w » OTF (y) 
Epa ~ peat (9) 
otw 04F(y) 
( oy4 | 2 Oy4 : i) 


L 


ee ee 
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es SO ee ee 


Ox? Oy? 4y? \a Oy? ys a Oy 
11 
Se R25 (on \2 a 
25 Gagrsal pyle 
The derivatives of F(y) are written as 
‘ Te uy 2uy 
yi =— = 
H h [cos h cos—> le ] 
gt \? cea, Be agt 
NES see — 
1 Gi) sin h + 2 sin h ‘fF 
ee a é ; cd T 2 h 
LP () [sin - — 8sin z i 
; I 2m y 3 3a Y 1 mY 
Py, == ||2 
2 h [ COS > cos— > cos > \ 
Seg eae age TC ea 3ay Le eat. 
Fy = (j) [ 4 sin eg Silent Fail l (12) 
eee a\t : y 81 3uy ul Tay, 
Ebi ae (F) [16 sin 2 sin h 2 sin i |: 
eS > ae ny 2a 
i = h [3 cos— > 2 cos cos—> : 
seg ele gt Spat! a ty ‘ e] 
F = (3) | 9 sin h + 8 sin as pee SiS , 


4 E eet 
Ry = (5) [st sin 2% — 128 sin 5 8 sin 7]. | 


If equations (12) are substituted into equations (7) through (11) and the resulting 
derivatives are substituted into equation (1), there results 


A ou + Ba,+ Ca, = 0, at y= 


a8 


2 
2h 
Ju+E a+ Fa,= 0, a i= 7, (13) 
3) 
Go+Hoa+T] «= 0, at y=, 
h\4 h\? : 
ie (x 4) + 72+ 22-00490 ( + 09101 32.6% 
h\2 
= 394-7860 — 21-5740 G 4 ; : (14) 
Nee h\2 | 
C =-—-x7? ( 2) — 799-4416 — 179-9193 (x | 0-101 32 6; 
a 


~ 
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2 
=) + 44-41342 + 5: 26442 G *) — 0:07599 B, ] 


J = 234221 (i - 


h\2 
Eb= 234011 (x -\ = $1-4246De 4113184 (x 4 + 0-07599 B, 
a 
h\2 
F = —592:1790 + 2178418 (i 4 
(14) 
h\4 h\2 
G = 2.35332 (: + + 85-93610 + 1-77300 (: <) — 0-12230 B, 
h\4 h\2 
H = —3-32809 ( 4 — 444-04924 — 60-45476 (i 3) + 0:17296 B, 


Gea igs h\2 
2S Soe (i *) + 644:24798 + 123:31292 (x 3) — 0-:12230 8 | 


and where eae 
Bln ey 5 (15) 
Values of f for various ratios of h/a may be determined from the condition 
pe tae canal 
[| eee pea 2 0\r (16) 
Ga tiga 


Figure 1 shows in graphical form, on semilog paper, the relationship between y 
and h/a for a practical range of h/a values, y = 61/2. For the particular case of a 
right isosceles triangle, from equation 16, y = 24-028. The second mode of vibra- 
tion of a square plate corresponds to the fundamental mode of a right isosceles 
triangle. The exact answer for the second mode of a square plate is y = 24-6741), 
Thus it is seen that the above solution is in error by 2-61 percent in the case of a 
right isosceles plate. 


3. Clamped Plates 


Consider the plate shown in the lower part of Figure 1. The boundary con- 
ditions are 


Ow Ow 
ya rs (10) (ay ae: (35 Ee SS eprawe ope ee 
The deflection function 
W = (% Fy + % Fy) Fy , (18) 
where 
F, = 9? sin? =~ tesa sin” sin 3 ; z ‘ 
7 , 
F, = [1 (=: =) ] cosa SS 
CUE aly a OEP Oy 


1) S. TimosHENKo, Vibration Problems in Engineering, 2nd edition (D. Van Nostrand Co., 
New York., 1937), p. 423. 
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satisfies the boundary conditions. It is necessary to differentiate equation (18) 


Vee ee ee eae 


_ and substitute the proper derivatives into equation (1). Let 
4 Fly)=uyR+ ak. (19) 
3 Derivatives of w at * = 0 become 
3 (@),-9 = Fly), 
 ( 02w hy? 1 Hig ate 
Z Paes Fel p= 
cape Fy) (7) y? [2+ e! | 
Bio) os (ey sia | 
Vor*) 2-0 SN) aya feels (20) 
io _ 0?F(y) 
| Oy? )x=0 Te Oy?’ 
(57) ee Onna) 
Oy* }x=0 Pay Oye ’ 
04w _ [A\? 0?F (y) 1 iene OF(y) 4 kh? a? 
“Seago (7) [- Oy? oy? [2+ 4 |+ dy 7 [2+ 4 | 
i! ie taps 
+ (—6 Fv) =a [2+ 7 I}. 
_ The derivatives of F(y) are written as 
es Bes cn ae ne Ming! 
Fy ie Y* Sia + 2ysin an | 
eet, ur 200 ae Qt 47 - gt Y 
EF, 4 = y sin i —E 2 y? cos a= 2 Sint rae 
oe % 2a Fike ey ae, eae 2.90 -y 
gl! te 24 75 cos ae 32 Ga Y sin— 8 qe Y° COS — | 
ee my 3a | M4 UY Bizet 72] 
a [eos h cos—> | + 3 7 2G + 7; sin; ; 
E (21) 
es Quy ae ; See [ : Sark = | 
ee h [-sin 7 + 3 sin h -+- |cos os h 
y? 2 uY K 32 | 
- ( Dh? ) [eos fi 9 cos 73 , 
Raised OE sh Cag i PA si Ot ogi 
ey 72 cos > Sit rr a ae eeMies beat ead h 
“y2 m4 ay > 3 7 | 
ae (+ Vo ) [cos ; 81 cos h ; 


If equations (21) are substituted into equations (20) and the resulting derivatives 
are substituted into equation (1), there results 


tS 
> 


Pion Ota, 0", Dee, Ry S a= 0, at y = 3 43 (22) 


a 
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where 
h\4 h\2 
P= 118"435 (,) hk? T — 42-052 + 39-478 () W — 0-25000 6, 
a 
OF 736:35 0), ; 
h\4 h\2 t (23) _ 
R= 49-965 (;) hk? T + 627-583 + 19-742 () W — 0:-33333 B, / 
a! 


2 
S = —49-965 (=)’ k? T — 1130-014 — 93-761 (=) W + 0:33333 B | 


and where 


ht Case a oe ee eee 
B fidgies ht, Seren —- 
Values of 6 for various ratios of h/a may be obtained from the condition 
EG 
= 0. 
eS 


Figure 1 gives of y, where y = 61/2, for practical ranges of the parameter h/a. In 
the case of an equilateral triangle 


y= 72-41. (24) 


In order to check the above value, the lowest natural frequency for an equi- 
lateral triangular plate is determined by the finite difference method. An equi- 
lateral triangular mesh is used1). The difference equations that result from writing 
equation (1) in difference form may be written matrically as 


ia Th ie a (25) 
where 
A = matrix of coefficients of 7? V2 terms; 


Jeanie matrix: 


4¢D 
s = length of side of triangle; 
m = number of subdivision points along a side of the triangle. 


If one symmetric point is taken, » = 4, and the value of y is found to be 42:31. 
If eight symmetric points are taken, » = 10, and the value of y is found to be 
65.85. Interpelation for a better value gives 

A (10)? (4)* 
There is about 3 percent difference between the above value and the value given 
by equation (24). The true answer should be between 72.41 and 70-34. 


: 1) A ts JENSEN, Analyses of Skew Slabs, Eng. Exp. Sta. Bull. No. 332, University of Illinois, 
p. 53 (1941). 
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It is seen that the solutions presented in Figure 1 appear to be quite reasonable 
and the simple curve should enable the engineer to determine rapidly the lowest 
natural frequency for clamped or simply supported isosceles triangular plates. 


Zusammenfassung 


Es werden Naherungswerte fiir die tiefsten Eigenfrequenzen gleichseitig drei- 
eckiger Platten mit eingespannten und einfach gelagerten Randern berechnet. 


(Received: June 22, 1954.) 


Durchgang von Explosionswellen durch unterkiihlte Wolken 


Von RAYMUND SANGER!) und FRITZ SPRING?), Ziirich 


Die Verfechter des Einsatzes gew6hnlicher mit einem Sprengkopf versehener 


_ Raketen zur Bekampfung der Hagelbildung schreiben deren hagelverhindernde 
_ Wirkung zum Teil den vom Explosionsherd ausgehenden Schallwellen zu. Durch 


re eae 


die Erschiitterung, die die unterkiihlten Wassertrépfchen einer Wolke beim 
Durchgang der Schallwelle erleiden, soll diese das Gefrieren der unterkihlten 
Trépfchen veranlassen und dadurch die Vereisung der Wolke erzwingen. Die 


 Intensitaét einer von einer Explosion ausgehenden Explosionswelle nimmt aber, 
verstirkt durch die auftretende Absorption, mit der Entfernung rasch ab, so dass 


die Einleitung der Vereisung héchstens in der Nachbarschaft der Explosion allfallig 


- erwartet werden kénnte. Nath den heutigen Kenntnissen ist aber die Vereisung 
_ einer unterkiihlten Wolke in erster Linie durch Zahl und Art der vorhandenen 
 aktiven Gefrierkerne bedingt, so dass die der Explosionswelle zugeschriebene Wir- 


kung nicht schlechthin verstandlich ware. Da schliesslich nur das Experiment 
entscheiden kann, wurden einige abklarende Testversuche durchgefuhrt. 
In einer ersten Gruppe von Versuchen, die im Schiesskanal der Werkzeug- 


- maschinenfabrik Oerlikon Biihrle & Co. ausgefiihrt worden sind, ist in einer 


Kihltruhe durch Einhauchen von Atemluft eine unterkiihlte Wolke von — 15° er- 
zeugt worden und hierauf Sprengk6rper, bestehend aus einem Hexogon-Trotyl- 


gemisch (85/15) mit nur gasférmigen Verbrennungsprodukten, im Abstand von 
5m und erhoéht gegeniiber der Truhe zur Explosion gebracht worden. Bei allen 
drei Versuchen konnte man jeweils eine lebhafte Bewegung der Wolke beob- 
 achten, die sich aber rasch beruhigte, doch waren: nicht die geringsten Anzeichen 


einer Anreicherung der in einem Lichtstrahl gewohnlich nur vereinzelt aufleuch- 
tenden Eisteilchen wahrzunehmen. Erst nach 4—6 Minuten 'trat eine klar erkenn- 


- pare Zunahme der Eisteilchen ein; diese Erscheinung kann aber ihre Erklarung 


darin haben, dass die Abgrenzungen des Schiesskanals stark mit dem vom 


Kugelfang herriihrenden Quarzsand verseucht sind. Durch die Explosionswelle 


wird aber mindestens zum Teil der Quarzstaub in Bewegung versetzt, und nach 
Ablauf einiger Minuten kénnten tatsachlich Quarzteilchen bis zur Truhe gewan- 
dert sein, wobei zu bedenken ist, dass Quarzteilchen ahnlich wie AgJ schon bei 


1) ETH. 
2) Werkzeugmaschinenfabrik Oerlikon Bihrle & Cie. 


~ 
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verhaltnismassig hoher Temperatur sich gefrierkeimfahig zeigen. Bei zwei weiteren 
Versuchen wurde dem Sprengsatz noch Silberjodid in Pulverform beigegeben; 
zundchst verhielt sich die unterkiihlte Wolke wie bei den ersten Versuchen, doch 


trat nach etwa 10 Sekunden unversehens ein starkes Anwachsen der Anzahl der | 


Eisteilchen in Erscheinung, ein Zeichen dafiir, dass aktiv sich erweisende AgJ- 
Teilchen die Kiihltruhe erreicht haben miissen. Die geschilderten Ereignisse 
sind in einem Film festgehalten worden. 


Da das dem Quarzstaub zugeschriebene nachtragliche Auftreten vereinzelter _ 


Eisteilchen eine Unsicherheit in die Beweisfiihrung der Untersuchung bringt, 
sind in einer zweiten Gruppe die Versuche im Freien auf dem waldigen Vorge- 
lande des Gaslaboratoriums der EMPA wiederholt worden. Die unterkiihlte 
Wolke wurde zunachst gegen die infolge der allgemeinen Luftbewegung auftre- 
tende Randturbulenz durch Uberdecken der Kiihltruhe mit einem diinnen Tuch 
geschiitzt. Sofort nach der Explosion ist jeweils dieses durch eine Glasplatte, die 
die Beobachtung erlaubte, ersetzt worden. Die Beobachtung allfalliger Eisteil- 
chen erfolgte in einem Lichtstrahl unter etwa 45° Neigung zur Kiihltruhen6ffnung 
und mindesrens wahrend 10 Minuten. Die 4ussern Abmessungen der Truhe sind 
30X35X95 cm (Hohe). Wahrend der Messungen herrschte leichter Westwind bei 
einer Lufttemperatur von etwa 15° C (11. August 1954). 

Bei den ersten drei Versuchen erfolgte die Explosion in einem horizontalen 
Abstand von der Truhenmitte von rund 900 cm bzw. in einer Uberhéhung gegen- 
liber dem Truheneingang von etwa 330cm; beim vierten Versuch von 280 cm 
bzw. von etwa 390 cm. Die Temperatur der unterkiihlten Wolke war beim ersten, 
zweiten, dritten bzw. vierten Versuch —18, —19, —14 bzw. —19°. Die zugehorigen 
Angaben der Sprengkorper sind: 10 g Hexogen-Trotylgemisch; 10 g desselben; 
25 g eines tiblichen Hagelraketensprengsatzes (26,3 % Al, 73,7% KCIO,), bzw. 2g 
desselben. Trotzdem mindestens 10 Minuten beobachtet wurde, konnten keine Eis- 
teilchen tiber das gelegentliche Auftreten vereinzelter Partikel festgestellt werden. 
Die Messungen bestatigen damit die zu Anfang ausgesprochene Mutmassung, dass 
die von einer Explosion des Sprengsatzes einer Rakete ausgehenden Schallwellen 
unterkiihlte Wolken nicht zur Vereisung bringen, auch nicht in der Nachbarschaft 
des Explosionsherdes. 


Summary 


The passage of sound-waves through supercooled clouds was experimentally 
examined ; no influence on the supercooled water droplets could be detected. 


(Eingegangen: 15. November 1954.) 
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Grundlagen der Elektronenoptik. Von W. GLASER (Springer-Verlag, Wien 
1952), 699 S., 444 Abb., sFr. 124.-. 

Die Elektronenoptik ist ein junger Wissenszweig der Elektronenphysik, bei 
dessen Aufbau seit Anfang der dreissiger Jahre Prof. GLASER durch seine umfang- 
reichen theoretischen Arbeiten massgebend beteiligt war. So kénnen wir das 
heute vorliegende Werk Grundlagen der Elektronenoptik als eine Rekristallisation 
der seit iiber zwanzig Jahren geleisteten Arbeit betrachten. In der Tat konnte 
sich der Autor in der Darstellung seines Buches weitgehend auf eigene Original- 
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_ arbeiten stiitzen. Wir treffen auch hier eine klare und scharfe Prazision der Be- 
_ griffe und eine konsequente Durchfiihrung der Gedankengange. 
| Das Buch zerfallt in drei Hauptabschnitte. Der erste Abschnitt, der die Halfte 
des Buches umfasst, beschaftigt sich mit den Abbildungsfeldern und der Gauss- 
_ schen Dioptrik. Ausgehend von den auf die Elektronen in elektrischen und magne- 
tischen Feldern ausgeiibten Kraften werden die Elektronenbahnen und ihre 
_ Eigenschaften behandelt, wobei von Anfang an der Anschluss an die Methoden 
' der gewohnlichen Lichtoptik iiber das Prinzip der kleinsten Wirkung gesucht 
- wird. 
__ _Imzweiten Abschnitt des Buches wird der Boden der Gau8schen Dioptrik ver- 
_ lassen und die Abbildungsgesetze dritter Ordnung, die Seidelsche Dioptrik be- 
_handelt. Unter Heranziehung der Wirkungsfunktion (Eikonal) werden zuerst all- 
_ gemeine Integralausdriicke fiir die Bildfehlerkonstanten des rotationssymmetri- 
_schen Feldes abgeleitet, um dann einzelne Bildfehler eingehend zu diskutieren. 
Anschliessend werden zwei Kapitel der Berechnung der Kaustik gewidmet. Die 
beiden letzten Kapitel dieses Abschnittes beschaftigen sich mit den Ablenkfeh- 
lern und der Fokussierung der kreisf6rmigen Elektronenstrahlbiindel. 

Im dritten Abschnitt widmet sich der Autor der Behandlung der Strom- 
dichteverteilung bei der elektronenoptischen Abbildung, insbesondere der Be- 
handlung von Beugungs- und Interferenzerscheinungen, wozu man allein an Hand 
der wellenmechanischen Auffassung des Elektrons gelangt. Der Ausgangspunkt 
ist die Hamilton-Jacobische Differentialgleichung in einer verallgemeinerten 
Form, welche, mit dem Erhaltungssatz fiir Teilchenzahl kombiniert, sehr an- 

_schaulich zu der Schrédinger-Gleichung fiihrt. Durch die Integration der Schro- 
dinger-Gleichung fiir die paraxialen Strahlen wird die Méglichkeit einer objekt- 
treuen Abbildung auch wellenmechanisch nachgewiesen, um hierauf auch die 
Seidelsche Dioptrik wellenmechanisch zu begriinden. Abschliessend folgt noch 
eine eingehende Behandlung des Aufldsungsvermégens des Elektronenmikro- 

- skops. 

Es ist dem Autor glanzend gelungen, die Methoden der theoretischen Elek- 
tronenoptik ausfiihrlich darzustellen, was in dem Vorwort als ein erstes Ziel des 
Buches angegeben wird. Der zweite vom Autor ersehnte Zweck, dem Techniker 
ein Nachschlagewerk fiir die Dimensionierungsfragen der Elektronenstrahlgerate 

zu bieten, ist nach der Auffassung des Referenten mit geringerem Erfolg erfiillt. 

Fiir ein ernsthaftes Studium der theoretischen Elektronenoptik ist das Buch 
warmstens zu empfehlen. E. B. Bas 


Elektromagnetische Strahlungsfelder. Von H. Zunrt (Springer-Verlag, 
Berlin 1953). 473 S., 170 Abb., DM 53.50. 

Dieses Buch behandelt vor allem die fiir die Nachrichtentechnik interessanten 
Fragen der Theorie der elektromagnetischen Wellenausbreitung. Es beginnt mit 
den Grundgleichungen der Elektrodynamik, fiihrt tiber die Wellengleichung in 
verschiedenen Koordinatensystemen, optischen und energetischen Betrachtun- 
gen zu einem eingehenden Studium der heute bekannten Antennenarten. Diese 
werden von der Theorie her dargestellt und berechnet, die Berechnung jedoch 
so weit getrieben, dass sie auch fiir die Bediirfnisse des Praktikers unmittelbar 
nutzbringend ist. Die letzten Kapitel sind der Erfassung von allgemeineren Aus- 
breitungsbedingungen gewidmet (Erdkriimmung, Einfluss der Atmosphare, Hin- 
dernisse). Ein ausfiihrliches Literaturverzeichnis ist bis 1952 nachgefiihrt. Fur die 
konsequente Anwendung des Giorgischen MaBsystems diirften dem Verfasser 
speziell die Studierenden dankbar sein. 
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Das Buch hat vor allem den Vorzug, Antennen- und Ausbreitungsfragen der 
verschiedensten Frequenzbereiche einheitlich und theoretisch sauber zusammen- _ 
zufassen und gleichzeitig eine weitgehende Einfiihrung in die mit einzelnen Wel- 
lenlangen verkniipften Spezialdisziplinen zu bieten. Dem wissenschaftlich inter-_ 
essierten Hochfrequenzingenieur diirfte es bald unentbehrlich werden. T. Stutz — 


Tensor Analysis, Theory and Applications. By I. S. Soko_nikorF (John — 
Wiley & Sons, New York, 1951). 335 pp., 43 Abb.; $6.—. . 

Dieser in der von I. S. SoKOLNIKOFF herausgegebenen Sammlung «Applied | 
Mathematics Series» erschienene Band ist aus Vorlesungen hervorgegangen, die 
der Verfasser waihrend mehrerer Jahre an drei verschiedenen amerikanischen 
Universitaten gehalten hat. Das Buch kann Studenten der angewandten Mathe- — 
matik, Physikern und Ingenieuren als Lehrbuch zur Einfiihrung in die Theorie © 
und Anwendung der Tensoranalysis sehr empfohlen werden. Es zeichnet sich 
durch eine geschickte und iibersichtliche Darstellung aus; die Begriffe werden 
sorgfaltig eingefiihrt und die gewonnenen allgemeinen Ergebnisse auf einfache, 
instruktive Beispiele angewandt. Jeder Abschnitt enthalt Ubungsaufgaben und | 
gelegentlich Literaturhinweise fiir ein eingehenderes Studium. Auch die Aus-_ 
stattung des Buches ist vorztiglich. | 

Das erste Kapitel enthalt eine knappe Einfiihrung in die lineare und quadra- 
tische Algebra und Matrizenrechnung, wahrend das zweite sich mit der Tensor- 
algebra und der Tensoranalysis befasst (es wird die tibliche Schreibweise und nicht 
die Kern-Index-Methode verwendet) und bis zum Riemannschen Kriimmungs- 
tensor vordringt. Die vier folgenden Kapitel sind den Anwendungen gewidmet. 
Zunachst wird ziemlich ausfiihrlich die Flachentheorie (innere und dussere) im 
FE, behandelt und, soweit erforderlich, die Riemannsche Geometrie im F,. Als- 
dann folgt eine gute Darstellung der analytischen Mechanik, einschliesslich eini- 
ger Hinweise auf die Potentialtheorie. Im fiinften Kapitel geht der Verfasser 
kurz auf die relativistische Mechanik ein und leitet hier die Schwarzschildsche 
(statische kugelsymmetrische) Lésung der Feldgleichungen ab, um die Perihel- 
drehung einer Planetenbahn zu berechnen. Das letzte Kapitel vermittelt noch 
eine Idee von der Anwendung der Tensorrechnung in der Mechanik der Kon- 
tinua. Es wird der nichtlineare Fall dargestellt und sodann der Ubergang zur 
klassischen Elastizitatstheorie und Hydrodynamik vollzogen. E. Roth-Desmeules 


Die Laplace-Transformation und ihre Anwendung. Von P. Funk, 
H. Sacan und F. SELie (Deuticke, Wien, 1953), 106 S., 18 Fig.; sFr. 10.-. 

Es ist das Ziel der Verfasser gewesen, eine kurze Einfiihrung in die Theorie 
und Anwendungen der Laplace-Transformation zu schreiben, die eine Mittelstel- 
lung einnehmen soll zwischen: den. fiir Mathematiker und den fiir Ingenieure 
geschriebenen Biichern. Diese, vermittelnde Darstellung ist vorziiglich gelungen 
und sehr zu begriissen, da einerseits die mathematischen Monographien der La- 
place-Transformation fiir den Techniker oft zu schwer lesbar sind und anderer- 
seits die mehr auf die Anwendungen zielenden Lehrbiicher — haufig noch von den 
alten Heavisideschen Operatorenrechnung inspiriert — es oft an mathematischer 
Strenge fehlen lassen. 

Hervorgehoben seien die mathematisch einwandfreie Darstellung des Um- 
kehrintegrals und der Asymptotik der Laplace-Transformierten sowie die knap- 
pen und das Wesentliche treffenden Einfiihrungen in die Anwendungsgebiete 
(Elektrische Schwingungskreise, Warmeleitung, Regelvorginge, Stabilitat). 

E. Stiefel 
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_ Mathematical Aspects of the Quantum Theory of Fields. By K. O. 
FRIEDRICHS (Interscience Publishers, Inc., New York, 1953). 272 pp., $5.00. 

7 Der Autor sagt in der Einleitung, dass dieses aus fiinf Teilen bestehende, auf 
Vorlesungen im Seminar tiber mathematische Physik an der New York University 
basierte Buch «in erster Linie an Mathematiker gerichtet sei, die mit den Grund- 
begriffen der Quantentheorie eines einzelnen Teilchens vertraut sind und die 
Jernen wollen, welche mathematischen Begriffe in den einfachsten Problemen der 
‘Feld-Quantentheorie verwendet sind.» 

Es ist klar, dass das selbst ftir einen theoretischen Physiker, der Spezialist in 
der Feldquantisierung ist, schwierig lesbare Buch nur auf einen sehr kleinen Kreis 
von Lesern rechnen kann. Obwohl das Buch nicht beansprucht, einen Fortschritt 
in den grundsdtzlichen Fragen der Feldquantisierung zu bringen, muss dem Autor 
zugestanden werden, dass es hinsichtlich mathematischer Strenge eine gewisse 
Liicke in der bestehenden Literatur in niitzlicher Weise ausfiillt. Dies gilt insbe- 
sondere hinsichtlich der konsequenten Vermeidung der Einfiihrung eines Hohl- 
raumes von endlichem Volumen mit Hilfe einer direkteren Verwendung der be- 
kannten Darstellung der Feldoperatoren von V. Fock. Ausser der iiblichen «Teil- 
chendarstellung» der Feldoperatoren werden noch andere mathematische Dar- 
stellungsméglichkeiten fiir diese Operatoren angegeben, deren physikalische Be- 
deutung vielleicht noch zu priifen ware. 

Das Buch beschrankt sich auf den Fall linearer Feldgleichungen, wo noch 
keine wesentlichen Divergenzschwierigkeiten auftreten. Uber den kraftefreien 
Fall hinausgehend, werden hierbei insbesondere auch noch die Felder gegebener 
Quellenverteilungen behandelt mit Beriicksichtigung der Divergenzen der Sto- 
rungstheorie bei kleinen Frequenzen («Ultrarotkatastrophe»). W. Paul 


Physikalisches Worterbuch. Von W. H. WestrpHat (Springer-Verlag, 
Berlin 1952). 1628 S., 1595 Textfig.; DM 148.-. 

The book represents an encyclopedia in one volume covering the fundamen- 
tals of physics, together with related fields of mathematics, astrophysics, and 
engineering. It is a completely new work continuing the tradition of the Hand- 
wérterbuch by BERLINER and SCHEEL, the last edition of which appeared some 
twenty years ago. The book is mainly intended for physicists as a means of rapid 
information on subjects outside their own special field. There is no doubt that the 
W6rterbuch will satisfy a definite need in this respect. 

The bulk of the material is arranged alphabetically in encyclopedic fashion. 
In addition, a concise history of physics is given in an appendix, followed by a 
number of tables of units and conversion factors, dimensional relations between 
the various electromagnetic systems, atomic weights, etc. A brief account of the 
more important astronomical numbers is also included. Finally, a list of physicists 
is given. It is obviously very difficult to select such a list and even within the 
stated selection rules there are inconsistencies, for example, WENTZEL and WIc- 
NER are not included. Such a list of names and birth dates does not seem to be 
of any use and could well have been omitted altogether. 

Any review of the large amount of material in the book, consisting of some 
10,000 headings, will obviously be weighted by the interests of the reviewer. In 
general, the impression is very favorable indeed. The text is mostly clear and 
concise, the many illustrations excellent, and the information has been brought 
up to the date of the printing (1952) with the help of a short addendum. Thus, 
recent work on V-particles, the use of Cherenkov radiation for counting purposes, 
etc. is included. Experimental techniques are well described and technical details 


. 
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often given. The ever increasing field of, electronic measuring methods seems 
adequately covered. Kae a iS SS ee 

It is clear that a book of this scope cannot be complete, and that certain 
- omissions and some shortcomings in special fields will have to be tolerated. More 
serious omissions noticed by the reviewer are the following: The grand cano- 
nical ensemble and its use in fluctuation theory is not mentioned. The central 
limit theorem of probability theory, the generalized harmonic analysis of WIENER: 
et al. are missing, and the probabilistic information concepts of SzILARD and | 
SHANNON are not sufficiently covered. The discussions on the ergodic hypothesis, | 
on the H-theorem, and on determinism could be extended and deepened. The 
heading ‘Thermodynamics’ does not do justice to the subject; to judge the 
merits of kinetic theory, statistical mechanics, and phenomenological thermo-- 
dynamics in a general way is fruitless and very much time dependent. Compared — 
with other sections dealing with hydrodynamics, the section on ‘Turbulence’ 
seems too brief and incomplete. Very little is included on non-linear vibration and 
especially a general concept such as limit cycles is missed. In some cases it is felt — 
that the presentation could be improved by a combination of several subjects _ 
_ under a single heading: for example, the classification of differential equations | 
into elliptic, parabolic, and hyperbolic types is not given in the section ‘Partial - 
Differential Equations’, but separately. The concept of GREEN’s function appears — 
under the separate heading of ‘Sturm-Liouville Differential Equations’. It 
would seem more economical and direct to combine the contents of these separate | 
sections into an article on partial differential equations, with back references 
under the more specialized sections; this has been done in other subject matters | 
such as, for example, ‘ Ultrasonics’. : 

Compared to the scope of the work this and similar criticism which can be 
raised appears, however, quite minor and unavoidable in a first edition. Paper, 
printing, etc. conform with the known high standards of Springer. A wide distri- 
bution, which the volume most certainly deserves, is probably only hampered by 
the necessarily high price. H.W. Liepmann 


Leitfaden der Nomographie. Von W. MEYER zUR CAPELLEN (Springer- 
Verlag, Berlin, 1953), 178 S., sFr. 20.10. 

Im ersten Teil dieses Werkes werden die theoretischen Grundlagen der Nomo- 
graphie entwickelt, im zweiten sind sehr zahlreiche instruktive Beispiele und An- 
wendungen zusammengefasst und zum Teil bildlich dargestellt. Diese fiir ein 
Lehrmittel der praktischen Mathematik etwas eigenwillige Gliederung wird durch 
gelegentlich schon im theoretischen Teil eingestreute Beispiele aufgelockert. 

___ Es ist dem Verfasser gelungen, in knappen 92 Seiten einen recht umfassenden 
Uberblick iiber die in der Nomographie angewendeten Methoden zu geben, fiir 
den die Bezeichnung «Leitfaden» beinahe zu bescheiden klingt. Eingehender als 
iiblich werden auch Fragen der Genauigkeit erértert. Dem Praktiker wird die 
tabellarische Zusammenstellung zahlreicher Funktionstypen mit den zugehGrigen ° 
Nomogrammformen wertvoll sein. 

Zum vollen Verstandnis des Werkes sind elementare Kenntnisse der Diffe- 
rentialrechnung, der analytischen und projektiven Geometrie notwendig. Auch 
die gelegentlich knappe Darstellung setzt eine gewisse Vertrautheit mit mathe- 
matischen Gedankengingen voraus. Das neue Buch gehért deshalb in erster Linie 
in die Hand des Akademikers technischer Richtung, der es als Einfiihrung oder 
als Erganzung zu einer Vorlesung schatzen wird. E. Véllm 
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TRUB, TAUBER: ZURICH 


Theoretische Elektrotechnik 


in vier Banden. Von Prof. Dr. K. Kuhlmann, Ziirich, 


Bis jetzt sind lieferbar: Band Il 


Grundziige der Theorie des Wechselstromkreises 


und des einphasigen Transformators 
VIIl und 175 Seiten mit 121 Figuren. Ganzleinen Fr. 27.55, broschiert Fr. 23.40 (1947). 


Hauptabschnitte: Symbolische Darstellung stationarer Wechselstrome — Geometrie der 
komplexen Zahlen, symbolische Ortskurventheorie — Theorie des einphasigen Wechsel- 
stromtransformators. : 


«Der vorliegende Band wird dem Studierenden eine gute Unterstiitzung im Studium sein. 
Auch der in der Praxis stehende Ingenieur, der Freude an einer theoretischen Wissens- 
vertiefung hat, wird gerne zu diesem Werk greifen.» (Technische Rundschau.) 


Band III 


Grundziige der Theorie elektrischer Maschinen 
547 Seiten mit 328 Figuren. Ganzleinen Fr. 74.90, broschiert Fr. 68.65 (1951). 


Hauptabschnitte: Allgemeines — Die Wirkungsweise der elektrischen Maschinen — Die 
Wicklungstechnik der modernen Maschinen — Erzeugung magnetischer Felder — Induk- 
tivititsberechnung von Maschinenwicklungen — Induzierte elektromotorische Krafte — 
Magnetische Energie und Drehmoment elektrischer Maschinen — Verluste in elektrischen 
Maschinen durch Strom- und Eisenwarme — Das Rechnen mit symmetrischen Komponenten. 


«Wer in diesem Buch eine Theorie der elektrischen Maschinen sucht, wird sie finden. Uber 
das bescheidene Vorwort ,, Grundziige“ geht der Inhalt aber weit hinaus. Zwar sind die 
Hauptkapitel: Wirkungsweise der elektrischen Maschinen, Wicklungstechnik, Erzeugung: 
magnetischer Felder, Induktivitatsberechnungen, induzierte elektromotorische Krafte, mag- 
netische Energie und Drehmoment.elektrischer Maschinen, wirklich auf die Grundziige 
abgestellt. Sie enthalten aber so viele Einzelheiten, dass der an der Theorie der Maschinen 
interessierte Ingenieur der Praxis kaum ein Problem aus dem Gebiet des Elektromaschinen- 
baues finden wird, auf das er hier nicht eine Antwort findet... Das ausgezeichnet aus- 
gestattete Buch wird allen am Elektromaschinenbau interessierten Ingenieuren willkommen 
sein.» (Elektrotechnische Zeitschrift.) 


Zu beziehen durch jede Buchhandlung. Fordern Sie bitte Spezialprospekte an! 
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